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Einige Anwendungen funktionalanalytischer Methoden 
in der praktischen Analysis 


Zusammenfassender Bericht 


Von LotHar CoLiatz, Hamburg!) 


In der praktischen Analysis sind in neuerer Zeit die Bestrebungen, Fehler- 
abschatzungen fiir Naherungswerte und damit exakte Einschliessungen fiir ge- 
suchte Werte aufzustellen, bei weiten Klassen von Aufgaben erfolgreich ge- 
wesen. Insbesondere hat man bei nichtlinearen Aufgaben in vielen Fallen 
Existenz von Lésungen nachweisen und Schranken fiir die Lésungen gewinnen 
kénnen. Die Erweiterung des Wirkungsbereiches friiherer Ergebnisse beruht 
vielfach auf einer Anwendung funktionalanalytischer Methoden; allerdings 
bediirfen diese Methoden haufig noch einer Umgestaltung, und es bleibt oft viel 
Freiheit in der Art der Anwendung, so dass viel Erprobung notwendig ist, um 
die fiir numerische Zwecke jeweils giinstigste Art festzustellen. Daher ist es 
vielleicht verfriiht, jetzt schon einen zusammenfassenden Bericht tiber dieses 
Gebiet zu schreiben. Wenn hier trotzdem ein solcher Versuch unternommen 
wird, so ist der Gedanke bestimmend, dass durch diesen Versuch weitere Kreise 
durch die bereits erzielten Erfolge angeregt werden sollen, mehr als bisher 
Fehlerabschatzungen durchzufiihren, Erfahrungen zu sammeln und den Wir- 
kungsbereich der Methoden weiter zu vergréssern. 

Die Gedankengange werden manchem Leser vielleicht etwas abstrakt er- 
scheinen; deshalb sind Nr.1 und 2 absichtlich sehr breit geschrieben, zugleich 
auch, weil eine Fassung des allgemeinen Satzes tiber das Iterationsverfahren 
gegeben werden sollte, die unmittelbar die Anwendung auch auf nichtlineare 
Aufgaben gestattet. 

Vollstandigkeit beziiglich der Anwendungen Paeeatidnalyrischer Metho- 
den war nicht Ziel des Aufsatzes, insbesondere sind die Anwendungen, auf 
Eigenwertaufgaben, tiber die schon viel Literatur vorliegt (vgl. zum Beispiel 
[16], [22]?), fortgelassen worden. 


1. Einige Grundbegriffe aus der Funktionalanalysis 


Der im folgenden wiedergegebene, fiir das Weitere grundlegende Beweis ist 
im wesentlichen der Arbeit von WEISSINGER [25] entnommen, aber auf eine 


1) Forschungsstelle fiir praktische Mathematik der Universitat. 
2) Zahlen in eckigen Klammern beziehen sich auf das Literaturverzeichnis auf Seite 356. 
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Form gebracht, bei welcher die Voraussetzungen bei Anwendungen unmittel- 
bar nachpriifbar sind (vgl. die Zusammenfassung am Schluss der Nr. 2). 

Es sei ein (abstrakter) Raum R gegeben und f, /,, fy, ... Elemente von ihm. 
Die Elemente kénnen zum Beispiel Zahlen, Vektoren, Funktionen sein. Fiir 
je zwei Elemente /,, f, des Raumes (kurz fiir /,, f,€ R) sei ei «Abstand» 
| f: — fy|| als reelle nichtnegative Zahl definiert mit den Eigenschaften 
1. Symmetrie || /, — f2|| = |f2 —/A| 
2. Definitheit: In | /; — /,|| = 0 steht das Gleichheitszeichen genau fiir f; = /y. 
3. Dreiecksungleichung: 


l4z-fell S]A—fsl + fs — fel fir fy, fe fg€ R. (1) 


, 


Ist zum Beispiel R der Raum der im abgeschlossenen Intervall <a, b> mit 
a < } eindeutigen stetigen Funktionen /(x) einer reellen Veranderlichen x, so 
kann man als Abstand definieren: 


lf: — fell = Max pets f2(*) | ; 2) 


S736 (*) 


wobei W(x) eine in <a, b> festgewahlte positive stetige Funktion ist. 

Bei der Anwendung ist gewohnlich eine kommutative Addition von Ele- 
menten erklart und F beziiglich dieser Addition eine additive Gruppe; es gibt 
ein Nullelement 0 in R mit 0 + / = f fiir alle f€ R und zu jedem /f ein In- 
verses —/ mit f+ (—/f) = @. Dann kann man an Stelle des Abstandes etwas 
bequemer eine Norm definieren: 


7 =|7- 4]. (3) 


Die Norm ||/|| eines Elementes ist der Abstand vom Nullelement. 
Ein Teilraum F von R heisst vollstandig, wenn zu jeder Elementenfolge 
livdans: aus F mit 


,, im fm -< | = 0 


ein Grenzelement / aus F mit 
lim ||f — /, || =0 
n—> OO 


existiert. 

Bei dem obigen Beispiel der in <a, b> stetigen Funktionen ist bei dem 
Abstandsbegriff (2) die Gesamtheit der in <a, b> stetigen, den Ungleichun- 
gen (x) S f(x) S u(x) geniigenden Funktionen f(x) ein vollstandiger Teil- 
raum (vgl. BANACH[1]), wobei (x), #2(x) zwei festgewahlte stetige Funktionen 
mit (x) < u(x) sind. 
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Nun sei ein Operator T definiert, der den Elementen f von F eindeutig 
Elemente T/ zuordnet, die zu R gehéren. Der Operator T heisst lipschitz- 
beschrankt, wenn es eine Lipschitz-Konstante K gibt mit 


| ieee Ob eS Ae 2a hun allewiajs c i. (4) 


Bei dem Beispiel der in <a, b> stetigen Funktionen f(x), dem Abstand (2) 
und dem Operator T mit 


T (x) = f G(x, fl) a8, 6) 


wobei G(x, €) eine in a S x, € < b gegebene stetige beschrankte Funktion ist, 
gilt: 
eee eee ee hee | 
mit 
uy 
J 1G 91 WS 4g 


Waeece, W(x) : (6) 


2. Der allgemeine Satz tiber das Iterationsverfahren 


Es werde unter den Voraussetzungen der vorigen Nummer nach Lésungen u 
der Gleichung 


fd (7) 
gefragt und das folgende Iterationsverfahren angesetzt : 
dey OL ae (n= OFly Zr =(8) 


Es gehére u) zu F; man kann w,,, so lange bilden, als uv, in F bleibt; fir 
O<xn< m gilt 


|| nm — Un | = (aes pare Sb Un || SK | ¢ ma — Una ree | (9) 


<= K* un + — Mn+ firO Sr sn 
und nach (1) und (9): 


| m—1 | m—1 ps m—n ? 
| Um — Un a i Oe (Us14 ce, Me) = py | Usi4 — Us| — a, K | Un — Un | 
||S=” S=n t= 


Sie ( y x) || Hy — Mo] - 
Sl 
Nun sei K < 1; dann folgt 
Kk” 
|| uv, — u%|| . (10) 


{ 
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Es werde nun der Teilraum S (kurz als «Kugel» S bezeichnet) eingefiihrt, der 
alle Elemente 4 von R mit 


|*—-m|S i= 


K | Uy — Up || (11) 


enthalt. Wenn diese Kugel S im Teilraum F enthalten ist, so liegen nach (10) 
(fiir » =1) alle w,, in F, und die Iteration ist unbeschrankt durchfihrbar. 
Aus (10) folgt wegen K < 1 


he al tm — te = 0, 


also wegen der Vollstandigkeit von F die Existenz eines Grenzelementes u aus — 
F mit ' 
lim | «— u,|| =0. (12) 
Nach (1) und (4) folgt dann 
[Tu —wl s|Tu—Tu,| + Tu, — ul] =|Tu— Taq] + [to — 4| 


SK |u— ql + |l¢nsa — ull; 


fiir n > co geht nach (12) die rechte Seite gegen Null; es ist also || Tu — || =0 | 
und damit Tu=u; zugleich ist die Existenz einer Lésung « von (7) gezeigt. 
Es sei v eine etwa vorhandene weitere Lésung von (7), also Tv =v; dann 
gilt nach (4) 
|«—e| =|Ta—To]/sK|u—o], 


also wegen K <1 auch |v — v|| = 0 oder « = v; das heisst, die Gleichung (7) 
besitzt in F eine und nur eine Lésung. Diese liegt sogar in der Kugel S, denn 
nach (1) und (10) gilt 


|e — 45] Shea] + fem — Ma] S 4 — Hell + ae Mees — tol - 


Hier strebt || « — w,,|| nach (12) fiir m oo gegen Null, es ist also 
oe: 
|v — uy| tee Pall Uy — Uo| (13) 
das heisst, w gehért S an. 


Zusammenfassung 


In einem Raum R sei eine Gleichung (7) vorgelegt. Man definiere in R einen 
Abstand, der den drei in Nr.7 gestellten Forderungen geniigt, und wahle einen 
Teilraum F so aus, dass 
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1. F bet dem gewihlten Abstand vollstandig ist, 

2. Tf fiir alle f © F eindeutig erklart ist, 

3. Tin F einer Lipschitz-Bedingung (4) mit K <1 geniigt, 

4. mi dem gewiihlten uy auch u,= T uy und die ganze Kugel S nach (11) zu F 
gehort. 

Dann gilt: Die Gleichung (7) besitzt in F eine und nur eine Lésung u. Die Itera- 

tion (8) ist unbeschrankt ausfiihrbar, und die Folge u, konvergiert im Sinne von 


(12) gegen u. Gleichung (13) liefert eine Fehlerabschitzung fiir die Naherung uy, 
das heisst u liegt in der Kugel S. 


Zusatz: In manchen Fallen (zum Beispiel in Nr.5) kann man direkt er- 
kennen, dass alle u, in F bleiben; dann kann die Bedingung 4 fortfallen. 
3. Lineare und nichtlineare Gleichungssysteme 


Ein System von  Gleichungen fiir ” reelle oder komplexe Unbekannte 
Xa), -++> X(m) Sei auf die Form gebracht 


% (7) = PilH a), «++» X(n)) » @ =-1,...,.) “(14) 


wobei die g, gegebene, nach den %(,) in einem Bereich F des (x(y), ..-, %(n))- 
Raumes stetig partiell differenzierbare Funktionen seien, und es werde 


(15) 


gesetzt. Die bekanntesten Iterationsverfahren sind die Verfahren (vgl. von 
MIsES-GEIRINGER [14]): : 
1. Iteration in Gesamtschritten 


Xj) k+l Pi(% (a) bv veee Xn) x) : (k = 0,1,...) (16) 
2. Iteration in Einzelschritten 
in ger = Piha) wea --> KU-1) be “Ue ie) (=O Line) Le) 


Es bedeutet bei «(;), der erste, eingeklammerte Index die Nummer der be- 
treffenden Unbekannten und der zweite Index die Naherungsstufe bei der 
Iteration, bei der man von %()o, ---, %(n)o ausgeht; fasst man die %();, fiir 
721; ..., wz einer Spaltenmatrix (oder einem Vektor) x, zusammen, so hat 
man bei beiden Iterationsverfahren eine Zuordnung 


Naps = Tx, . (perl Last.) (18) 


332 LoTHAR COLLATZ ZAMP 


Nun lisst sich die Theorie von Nr. 2 anwenden, und es ordnen sich verschiedene 
bekannte, hinreichende Konvergenzkriterien durch geeignete Wahl einer Norm 
hier ein. 
A. Iterationsverfahren in Gesamtschritten 
1. Norm: Es werde fiir die Spaltenmatrix x mit den Elementen xj) die Norm 
| *|| = res | *a| (19) 


gewahlt. 
Fir zwei Spaltenmatrizen x und y mit den Elementen (;) bzw. yj) ist dann 
nach dem Taylorschen Satz 


|T«—Ty| = Mee lpi(%() — Psl¥~y)| S Man Dh Ay |X) — Yory| 
r=1 


= (ax pa a (Maz xm — Ym |) SK x — yl, 


wenn 
= Max )"a;, (20) 
1 r=1 


gesetzt wird. Die Bedingung K <1 ist hinreichend fiir die Konvergenz des 
Verfahrens und stellt das Zeilensummenkriterium dar. 
2.Norm. Es werde 


| = py xa (21) 


als Norm gewahlt. Jetzt ist 


|7x— Ty] = Slo X (ry) — Pi(Very Ns SS arly = "Vinyl 


j=l r= 


n n 
=D) 14m — Yin | DY Gir S Kx — y| 
r=1 j=1 


mit 
Kees Max 3" a,;,. (22) 
Tas, 


Die Bedingung K < 1 ist das Spaltensummenkriterium. 
3. Norm. Es werde 


lel=+ VS bol 23) 
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als Norm gewahlt. Unter Benutzung der Schwarzschen Ungleichung wird jetzt 


|Tae— Ty? = a | P4(%(n) — Ps(Vo) |? LL 137 Asp | Xn) — vol} 
<3 (3 at) (3 lx — sol?) =e 9S Sa 


1=1 7= 
Die Bedingung K < 1 mit 
K=+ >s ae, (24) 


ist das Erhard-Schmidtsche Kriterium. 
4. Norm. LoNsETH [13] schlagt noch andere Normen vor, zum Beispiel 


bln 
| *|| = &, mo? mite pe, ls (25) 
mm 


oder MORGENSTERN [15] noch allgemeiner 


Pp 7 ae q 
| *|| = P? Peas S |xol* 
s 


j=v+l1 


mit f>1,¢g>1,1S sn, P>0,Q->0. Die Lipschitz-Konstanten lassen 
sich bei der Norm (25) mit Hilfe der Hélderschen Ungleichung 


Sendai s |/ Se tovle f/ Sse 


mit (pf — 1) (g — 1) = 1 berechnen. (25) enthalt als Spezialfalle bzw. Grenz- 
falle (21) fir p > 1, (23) fiir = 2 und (19) fiir > o~. 


B. Iterationsverfahren in Evnzelschritten 


Es wird jetzt die Vorschrift (17) und (18) verwendet, und als Norm werde 
(19) benutzt. Wendet man den Operator (18) auf zwei Spaltenmatrizen ~, y 
mit den Komponenten %(;) bzw. (;) an und setzt man 


lay —ygl =i. Tay —-—(T9)ml=a,-D= ee d; (26) 


(bei Tx bzw. Ty bedeutet der eingeklammerte Index 7 wieder die 7-te Kom- 
ponente), so ist 


os Say vot Said oe (27) 
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Nun fiihrt man (nach SASSENFELD [19]) die Gréssen ein 
n j-1 ‘ n ‘ 
Bi= D7 are; Bs= D 452Bet Df Ma (fir j = 2,.... 0); K=MaxB;. (28) 
s=1 s=1 s=]j 


Dann beweist man e; < D f; leicht durch vollstandige Induktion: Fir 7 - 1 
ist sie nach (27) richtig; nun sei sie bereits fiir ¢= 1p2, .«, &— lachtig, 
dann hat man fiir 7 = & nach (27) 


k-1 n 
ty = D/ tts D Bs + D teeD—D fr. 
s=1 s= 


wie behauptet. 
Mit der Norm (19) ist also 


|Tx—Ty| =Maxe;SKD=K||x-y]|, 
] 
das heisst, das nach (28) berechnete K ist als Lipschitz-Konstante verwendbar. 
Wenn das Zeilensummenkriterium erfiillt ist und also nach (20) 


Max 7 aj,= <1 (29) 
J r=1 


gilt, kann man etwas grober (dafiir aber bequemer, da man die Berechnung 
der f; spart) abschatzen, indem dann 


j-1 n 
Bsus § SD) 44+ > 4,548 
sS=] s=] 


durch Induktion beweisbar ist und damit einfach « als Lipschitz-Konstante K 
benutzt werden kann. 

WEISSINGER ([25], S. 199-201) bewies die Konvergenz des Iterationsver- 
fahrens in Einzelschritten und die Giiltigkeit der Fehlerabschatzungen auch 
fiir die Falle, in denen die in (22) und (24) gegebenen Gréssen K < 1 sind. 

Fiir den vielleicht wichtigsten Fall des Kriteriums K <1 nach (20) gilt die 


Zusammenfassung 


Beim nichtlinearen Gleichungss ystem (14) werde ein Bereich F im (x1), ..+, X(n))- 
Raum gewihlt und die Grissen aj, nach (15) ermittelt. Fiir das Gesamtschritt- 
verfahren (16) berechnet man K stets nach (20), fiir das Einzelschrittverfahren 
kann man zwar K ebenfalls nach (20) bilden, genauer aber verwendet man dann (28). 
Ist K <1, so ist die Kugel S mit xi; als laufenden Koordinaten gegeben durch: 


K 
Mes [4 — *eu| S Ti whee 1%) nea — Xap al - (30) 
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Gehort diese Kugel S ganz zu dem gewdhlten Gebiet F, so konvergiert beim Itera- 
tuonsverfahren x,;), fiir k oo gegen das einzige in F vorhandene Losungssystem 
von (14), und dieses liegt sogar in S, so dass (30) zugleich eine F ehlerabschatzung 
darstellt. 


Spezialfall des linearen Gleichungssystems 
Lautet (14) a 
DES (iieen Lace fh) 
k=l 
mit konstantem «;,, so braucht man den Bereich F nicht zu wiahlen; die Vor- 
aussetzung, dass die Kugel S zu F gehort, fallt fort, da in (15) 


5 I 
O59 


konstant. werden, wobei angenommen ist, dass man das Gleichungssystem 
nach den Hauptdiagonalgliedern auflést, um es in die Form (14) zu bringen. 
Es ergeben sich dann unmittelbar die verschiedenen Kriterien in der gelau- 
figen Gestalt. 
Zwei Beispiele sollen die bequeme Anwendbarkeit der Formeln zeigen. 
Beispiele 
I. Gleichung 1 + z = e*; Iterationsvorschrift : 


Sea = Gee) me, (Zz) "In (1 + 2). (k= 0, 1, ...) 


Einige Iterationen ergeben: 


24 In (1 + 2%) 
2ni 1,85 +7,704 eon in = dela) 
1,85+7,70 1 2,105 +7,499 i Ber Aeer (hes. 0 os I 
2,10 +7,50 i 2,0944-7,462 1 To natn —0, 1,6.) 
gy = 2,09+7,46i | 2 = 2,08873+7,46128 1+ 200i (n =0, +1,...) 


Als Bereich F werde Rez > 0,Imz = 2a und dort als Wert des Logarithmus 
weiterhin der Zweig gewahlt, der sich stetig an den Wert 


In(2%7) =In2a+i> 
anschliesst. Dort ist 


a rit eee (4 eld OST 
jt+e2, = 
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als Norm werde der gewodhnliche Betrag || z|| = |z| gewahlt; man hat nur 
noch nachzupriifen, dass die «Kugel» S, hier der Kreis |z—2,| <7 mit 
y =[K/(1 — K)] |z3 — 24| = 0,00034 dem Bereich F angehort. Das ist der Fall, 
also gibt es (bei dem gewahlten Zweig des Logarithmus!) genau eine Losung z 
der Gleichung in F, die sogar in dem Kreis | z — z4| < 0,00034 liegt. 

II. Die beiden Gleichungen fiir x, y: 


x?+y?—-2=0, 3xy—yi—1=0 


werden umgeschrieben: 


2 


1 pea 
x = 93(%, ¥) - 5 a4? y = Pa(x,; y= V2 aaa 


Nach kurzem Iterieren kommt man etwa zu den Werten: 
%, = 0,755; y,=1,195;- —,(xz,'¥,)'= 0,75495, oa, ¥5)— 11 ese 
Man wahlt daher zum Beispiel als Gebiet F 


13195, S-ves Oss tO ote oes Oto 


dort ist 
ay, = Max ae = 0; ayy 0322 = Bs 
inF | 0% 
4 K 


Die maximale Anderung @ betragt 0,00085. Die «Kugel» S lautet nach (30) 


|x — 075495] < =2~ = 0,041; | y — 1,19585| < 2% = 0,00041 . 


De 


Da S ganz zu F gehdort, stellen diese Ungleichungen zugleich eine Fehler- — 
abschatzung fiir die einzige in F vorhandene Lésung x, y des Gleichungs- 
systems dar. 


4. Das gewohnliche und vereinfachte Newtonsche Verfahren 
bei nichtlinearen Gleichungssystemen 


Es sei das Gleichungssystem fiir 7 (reelle oder komplexe) Unbekannte 
X(y)) +++) X(m) Vorgelegt 


Fil Xajn <0, Mi) =O , (f <1 2y ceae) 


wobei die /; in einem Bereich F des x;,)-Raumes gegebene, nach allen Argu- 
menten stetig differenzierbare Funktionen seien. Beim gewéhnlichen Newton- 
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schen Verfahren berechnet man aus der k-ten Naherung die (k + 1)-te Nahe- 
Tung %(;) x4, durch Aufldsen des linearen Gleichungssystems 


+ Se ie) ou Te yveene Ie OS Ie a k= OS 1, ae noe) 


oder ktirzer in Matrizenschreibweise 


G(x) (%n+a — Xx) = — fe (33) 
mit 
0 i) 
mye fale) ae ro 
Xx == ) ie == 0.00.0 | 5 G | onelueinejielle) (elise) Wiel ele) sits . (34) 
Ofn Ole 
X (nye in(%x) 0% (4) O%(n) || 


WILLERS [26] weist darauf hin, dass es fiir die Rechnung bequemer ist (und oft 
ebenso rasch zur gewiinschten Genauigkeit fiihrt), in (33) nicht bei jedem 
Iterationsschritt G(x;,) neu zu berechnen, sondern mit einer festen Matrix G(z), 
etwa G(x), zu arbeiten. Zudem kann man dann die reziproke Matrix G—! ein 
fiir alle Male berechnen. Dieses Verfahren werde vereinfachtes Newtonsches 
Verfahren genannt. (Vgl. [6], dort auch eine Fehlerabschatzung fiir Nullstellen- 
bestimmung bei algebraischen Gleichungen mit Hilfe des Hornerschen Schemas.) 
Die folgende Berechnung der Lipschitz-Konstanten und die damit gegebene 
Fehlerabschatzung gilt auch fiir das gewohnliche Newtonsche Verfahren, da 
man stets den letzten Schritt des gew6hnlichen Newtonschen Verfahrens als 
Anfangsschritt des vereinfachten Newtonschen Verfahrens auffassen kann. 


Es ist dann 
Kya = %, — G*(z) fp = T %, . (35) 


Hierdurch ist eine Transformation T erklart. Sind x,, x, zwei verschiedene 
Wertesysteme, so ist 

Pipl te [Ge (2) G(Ex)] (%% — %x) » 
wobei E die n-reihige Einheitsmatrix und &, gewisse Zwischenstellen zwischen 


x, und x, bedeuten. Man kann nun wieder verschiedene Normbegriffe wie in 
Nr.3 verwenden. Mit der einfachsten Norm (19) wird 


WT Xnaa— T Xess | a0 | ee ag 


wenn 
k= Max) Ajn (36) 
J r=1 


ZAMP IV/22 
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gesetzt wird; dabei ist A = E — G~1(z) G(g,) eine Matrix mit den Elementen 
a;,, und A;, sind obere Schranken fiir die Betrage von a;,. Es ist G~*(z) eine 
Matrix mit festen Zahlen, aber in G(&,,) hat man die &, den Bereich F durch- 
laufen zu lassen. Bei diesem vereinfachten Newtonschen Verfahren braucht man 
zur Abschatzung nur die ersten partiellen Ableitungen der /;, wahrend man 
beim gewohnlichen Newtonschen Verfahren zur Abschatzung die zweiten par- 
tiellen Ableitungen mit heranziehen muss (vgl. KANTOROVITCH [11], S. 88). 


5. Anfangswertaufgaben bei gewohnlichen Differentialgleichungen 
Die Differentialgleichung 
y’ = f(x, 9) | (37) 
mit dem Anfangswert y(%») = V9 werde als Integralgleichung 
v(x) = yo+ | fle, ve] at (38) 
geschrieben und als Spezialfall der Volterraschen Integralgleichung 
y=Ty mit Ty(s)= y+ | Gls,t, yO] a (39) 
aufgefasst. Der Kern G(s, ¢, vy) sei im Bereich H: 
|s—S9| Sa*, |é—59| 5 a*,. |y—%| Sb 
integrierbar, beschrankt durch |G| < M, und geniige dort einer Lipschitz- 
Bedingung 


G(s, t, i) — G(s, t, V9) 


My ae 


CS ba (40) 


Es werde noch a = Min(a*, b/M) gesetzt. Zur Anwendung des Satzes von Nr.2 
wird der Teilraum F als die Gesamtheit der in |s — sy| <a (Intervall /) 
stetigen Funktionen y(s) mit | y(s) — yp| < b gewahlt, der nach Nr.1 bei dem 
Abstand (2) vollstandig ist. Geht man bei der Iteration 


Vn+1(S) =T ¥n(S) (n = 0,1, asl 
von einer Funktion yo(s) mit ¥o(so) = Vo und | yo(s) — ¥o| Sb aus, so bleiben 


alle y,(s) in F (es tritt der Zusatz in Nr. 2 in Wirkung). Nun ist noch zu 
priifen, ob in der Lipschitz-Beschrankung (4) des Operators K < 1 ausfallt. 
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Es wird, wenn y(s) und z(s) zwei Funktionen aus F mit y(s9) = 2(so) = Vo sind: 
(Lee 
W(s) | 


1 
S Max weg | 1Gls. 4, 9) — Gls, 4,2)| a 


| Ty — Tz| = Max 
if 


s (41) 
L 
= Say cz 
= ees wie | | z| dt 
E AY 
< Max wy | ¥O at|y—2|=Kly—all | 
mit ; 
[we dt 
K=L Vee 2 Ws) (42) 


MorGENSTERN [15] verwendet als W(s) die Funktion e** und zeigt, dass man 
fiir s > sy mit A > L stets K < 1 erreichen kann; denn es wird 

Ss 
en ea eae CaS 2 


A 


~| 


Man erkennt so zugleich die eindeutige Losbarkeit von (37) bzw. (39) am Inter- 
vall J unter den getroffenen Voraussetzungen; fiir numerische Zwecke ist es aber 


oft giinstiger, andere Funktionen W(s) zu verwenden. 
Beispiel: Nur zur Erlauterung sei ein ganz grobes Beispiel genannt. Bei 


y=ix => (8+), y(0) = —1 


werde der rohe Ansatz y)(x) =—1+a«4x gemacht und « so gewahlt, dass 

Vy, — Vo = 0(x) klein ausfallt. Hier ist 

14+ o? 
fe) 


x3 


Vay 1+ 5-3 + 


und es werde « = 0,45 gesetzt. Als Teilraum F wahlen wir die in 0% 51 
stetigen Funktionen y(x) mit —1,02 Sy 30; dort ist |df/dy| = |v| SL=1,02. 
Fir W(x) = 1+ 4 x? ist 

| W(t) dt 

0 val 7 


<0,1> W (4) eis, 
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also wird nach (42) 


0,254 
3] 


K =~ 1,02 = 0,476; weiter ist || y, — yo|| = 


, 


damit nach (13) 


0, na K 
|y— nl s-3™*- 7 = 0.00462, 


mithin hat man als « Kugel» S 
ly — yo] S 0,00462 (1+42%*%) fir OS x1. 


Da S in F liegt, stellt diese Formel zugleich eine Fehlerabschatzung dar. 

Die Ubertragung auf Systeme und damit auf Differentialgleichungen héhe- 
rer Ordnung ist leicht und zeigt eine grosse Analogie zu den Gleichungs- 
systemen in Nr.3. Die Funktionen yi), ..., ¥(,) werden zu einer Matrix y zu- 
sammengefasst, ebenso die Funktionen /,, ..., /, zu einer Matrix f und die 
Anfangswerte (;)(%9) = (4) 9 Zu einer Matrix yo: 


= Ie, V1)» eeey Y(n)) (43) 
oder 


x 


y=Ty mit Ty=y 9+ / {[E, My (t)] at . (44) 


Hiermit ist zugleich ein Operator T festgelegt. Die Funktionen f; mégen in 
einem beziiglich der yj.) konvexen Bereich H des (x, Y(1)> +++» V(m))-Raumes 
Lipschitz-Bedingungen pba ee 


| fj(x, Vajr eres Yin) — Tal, Z(q)y ++, 2 n)) | SJ) Lals) (x) — 2a) | « (45) 
7. Norm: 
i) 
yl] = Mae (Max Iv {i) merely. (46) 


Dabei ist J ein den Punkt x9 enthaltendes Intervall der x-Achse und W; (x) 
sind (fir 7=1, ..., ) festgewahlte, in J positive stetige Funktionen. Fiir 
zwei Matrizen y, z mit y(x9) = 2(x%») = Vo gilt dann 


fey Ein) [Yay — %a| 4 


1 


|T y—T 2 | S Max Max *** eee <|y-—2|K 
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mit 


K = Max Max 
i j W;(¥) 


(47) 


Nimmt man die L;, als Konstante und alle W;(x) einander gleich: W,(x) =W(zx), 
so wird 


j W(t) dt 


K = @ Max L;, mit @®= Max 
i Pa é pie 


(48) 


(*) 


Bei K hat @ als Faktor die maximale Zeilensumme in der Matrix der Lipschitz- 
Konstanten. 
2. Norm: 


ie (x) 
|] = Max (a4 >» aes oo) 
Jetzt wird entsprechend 


[Bie |%e) — 20y| @ 


|Ty—Tz| Ss eae ° Wi) 


Setzt man die W;(x) = W(x) einander gleich und fiihrt die maximale Spalten- 
summe in der Matrix der Lipschitz-Konstanten ein: 


= Max Y Lal) , 

2 ee 

so wird 
so | 519%) ~ “|| uy 
fos ah olan %o aE <= K|ly—2| 

mit 

[so W(t) at 

K= Max =" (49) 


i W(x) 


Wieder kann man bei beiden Normen W(x) = e** mit geniigend grossem A 
setzen und damit in einem hinreichend kleinen Intervall Existenz eines L6- 
sungssystems und Konvergenz des Iterationsverfahrens sicherstellen; fiir nu- 
merische Zwecke kommen oft auch andere Funktionen W;(x) in Frage. 
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Es interessiert noch die Anwendung auf das System 
Yin = Vos (ftir 7 = 1, 2,...,%- 1), omy = 8(% Yay +++» Vem)» 
das heisst auf die Differentialgleichung m-ter Ordnung 
yin) = g(x,y, y's... YO"). 
g geniigt hier einer Lipschitz-Bedingung: 
| (3, Sint, OR) ae BU ape gia-)) | Ly |v") = 22-9], (50) 


wahrend die iibrigen 


1 firk=7+1, (7=1,.:.,”#—1) 
Vhs. = 
* 0 sonst, das heisst firk +7+1, (7 =1,...,2—1) 


sind. Beide Normen sind anwendbar. 


6. Randwertaufgaben 


Fiir eine Funktion der unabhangigen Veranderlichen x,, ..., %,, sei in 
einem Bereich B des (x, ..., %»)-Raumes eine Differentialgleichung 
L[u] = 9(%4,---. %m,¢) in B ! (51) 


und lineare Randbedingungen 
U[u]=y,(% +.) %m) auf He (uw = 1,..., &). (52) 


vorgelegt. Dabei seien L{w] und U,{w] lineare homogene Differentialausdriicke _ 
in # und den partiellen Ableitungen von « mit Koeffizienten, die von x,, ..., Xm 
abhangen kénnen, und @ eine gegebene, etwa stetige, nach wu stetig partiell 
differenzierbare Funktion ihrer Argumente. J’, sind gegebene (m—1)-dimen- 
sionale Hyperflachen, gewohnlich Randflachen von B. 


Nun werde die Annahme getroffen, dass die lineare Randwertaufgabe fiir 
eine Funktion w 


LI) ee Fi tert a en cats 
(53) 
Ulw]=0 auf I',(u = 1,..., 2) 


bei beliebig vorgegebener stetiger Funktion r(x, ..., *,,) stets eindeutig lésbar 
sel. Das ist zum Beispiel der Fall, wenn eine Greensche Funktion G(x;, &;) 
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existiert, welche die Randwertaufgabe (53) in der Form 
w(x) = [ Gls, &) IE, a8, (54) 
B 


lést. Dann ist auch die Randwertaufgabe L[w] =, U,[w] = y, stets eindeutig 
lésbar und das Iterationsverfahren 
Lina] = gle, ty) in B | 


(55) 
OU [4nia] = 7, auld (i= I ky 0, 12). 4, | 


bei der man von einer gewahlten Funktion u,(x;) ausgeht, unbeschrankt durch- 
fihrbar, vorausgesetzt, dass p(x;, w) fiir alle w = u,, definiert ist. Fiir die Uber- 
legungen (nicht fiir die numerische Rechnung) werden die Randbedingungen 
homogen gemacht. Es sei w eine spezielle, die Randbedingungen (52) erfiillende 
Funktion, dann setzt man 


uU=U+N, tn=Ut Ny, (56) 
und definiert zu einer Funktion f/ eine Funktion g und damit einen Operator T 
mut = J}: 
Lg] = v(x; u+f/)—L[u] in B 
U[g]=0 at I (aL, a ee) 


(57) 


Dann gilt 
Mou=1n, N=ITn, *—%,=4—-%, (n=90,1,...). (58) 


Die Fehler von wu, und 7, sind also dieselben, man kann daher die Aufgabe 


fiir 7 weiter behandeln. acs 
Um den Satz von Nr. 2 anwenden zu kénnen, wird die Lipschitz-Konstante 


K von (4) gebraucht. Seien /,, /, zwei Funktionen mit g;= Tf; (¢ = 1, 2), 
so wird 


Igy ad 82] = 9(%;, u+ hi) — P(%s, u+ fe) , U, [1 TF | =a Up. (59) 
Nun wird ein Bereich H des (x1, .-., %n, “)-Raumes gewahlt, in welchem 
| o(%5, 4) — 9(%s, u*)| = N(x;) |u — u* | (60) 


fiir alle x;, w, w* in H gilt. (Da @ als nach wu stetig differenzierbar vorausgesetzt 
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war, kann man N = Max |0q/du| setzen.) Aus (54), (59) und (60) folgt daher 


gi — b= [Glen &) [Es + hy) — wes + fa] Os, (61) 
B 
len — eo] Sf NE) [Gln )| | — fal a8 (62) 
B 
Wahlt man als Norm einer Funktion 
1 #| If 
|f|| = Obere Grenze Wx)” (63) 


wobei W(x;) eine in B festgewahlte stetige positive (eventuell auch nicht- 
negative) Funktion ist, so folgt 


lex — gel =KlA-fl (64) 


mit der Lipschitz-Konstante 
| N(E) |G(x;, &) | W(E) 48; 


B 
K= Obere. Grenze Wie) : (65) 


Es sei auf einen Fall hingewiesen, bei dem man ohne Kenntnis der expliziten 
Form der Greenschen Funktion zu einer schnellen, aber meist nicht so genauen 
Ermittlung der Lipschitz-Konstanten kommen kann. Wechselt die Greensche 
Funktion G in B nicht das Vorzeichen, ist etwa G = 0, und gibt es eine eben- 
falls in B das Vorzeichen nicht wechselnde Eigenfunktion z(x;) (etwa z; = 0) 
von 


[[z]=Az inB, Uf[z]=0 auf I, (66) 
zum Eigenwert A = /,, so kann man W/(x;) = 2(x;) wahlen, und nach (54) gilt 
[Gln &) 2G) dy = Fal); 

es folgt einfach [ 


K = 2 maw, (67) 


(Verschiedene Zahlenbeispiele bei [5], Z. angew. Math. Mech. 33, 116-127 [1953], 
und ein Beispiel hier in Nr. 11.) 
7. Allgemeinere nichtlineare gewoéhnliche Differentialgleichungen 


Ist die Differentialgleichung (51) nicht nur in «, sondern auch in Ablei- 
tungen von w nichtlinear, so kann man manchmal ahnlich wie in Nr.6 vor- 
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gehen. Das sei kurz an dem Beispiel der gewéhnlichen Differentialgleichung 
k-ter Ordnung 

L[u(x)] = o(x, u, uv’, ..., u) (68) 


und Randbedingungen an zwei Stellen x = a und x = b (b > a, s < k) 


U4] =v, (u=1,..., ) 
genannt. 
L{w] sei ein linearer homogener Differentialausdruck in u, die Randbedin- 
gungen seien linear. m gentige einer Lipschitz-Bedingung: 


| p(x, “, uy, Gor, u\)) a p(%, Up, Us, seey uy) | 
: (69) 
= Dd? A(x) ul — uf)| mit A(x) 20 
a=0 


Wieder sei die Randwertaufgabe (53) durch (54) lésbar und (54) sei s-mal nach x 
differenzierbar : 


b 
w(x) = [ G\x, B)r(@) dé mit Go SEs), (70) 
definiert man einen Operator T mit g = Tf durch 
Hgl= 9% tf. f), Ulglare Wad. k), 
so erhalt man bei der Norm 
1 < . 
|| = hee Geis Way Dy Are TP)! (71) 
in (a, a0 


nach kurzer Rechnung (genauere Durchfiihrung bei [5], Z. angew. Math. Mech. 
33, 124 [1953]) fiir die Lipschitz-Konstante K von T 


KS Obs rene en Ty A fia x, &)| W(é) d (72) 


8. Operatoren monotoner Art 


Der Raum R sei nun ein reeller Raum, das heisst, seine Elemente seien 
reelle Gréssen, etwa reelle Zahlen oder Matrizen mit reellen Elementen oder 
reellwertige Funktionen reeller Veranderlicher. Der Raum sei halbgeordnet (vgl. 
zum Beispiel KANTOROVITCH [11]), das heisst, Beziehungen / Sg oder f/<g 
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seien fiir gewisse Paare von Elementen erklart und sollen die gewohnliche Be- 
deutung haben. Wieder sei ein Operator T gegeben, der den Elementen g eines 
Teilraumes F eindeutig Elemente T g in R zuordnet. Es sei nun f ein gegebenes 
Element in R, und es werde nach Lésungen u von Tu =f gefragt. Diese Auf- 
gabe heisst «von monotoner Art», wenn aus Tv S Tw fiir beliebige Elemente 
v,w aus F folgt v <w, und T heisst dann «ein Operator monotoner Art» 
(genauer: die Aufgabe ist als Aufgabe «von monotoner Art» geschrieben). 

Ist eine Lisung u vorhanden, so ergibt sich die Moéglichkeit einer Einschlies- 
sung. Sind v, und vg zwet Naherungen mit 


Tite ef ae Ete (73) 
so folgt 
Uy Sts S& Ue. (74) 


Man kann dann zum Beispiel nach Art der Relaxation versuchen, von einer 
Naherung v ausgehend, den «Defekt» d(v) = Tv — f durch Anbringen kleiner 
Korrekturen an v dem Nullelement zu nahern; gelingt es, dabei d = 0 zu 
erreichen, so weiss man v = wu, und entsprechend gilt v S uw bei d SO. Not- 
wendig ist aber bei dieser Methode die Kenntnis von der Existenz einer L6- 
sung #; man hat dafiir drei Wege zur Verfiigung: 

1. Oft gelingt es, mit Hilfe des Iterationsverfahrens und eines eventuellen 
ganz groben Iterationsschrittes nach dem Satz von Nr. 2 die Existenz einer 
Lésung zu zeigen (vgl. hierzu auch Nr.12 und Nr. 11, Beispiel I). 

2. Die Erweiterung des Browerschen Fixpunktsatzes fiihrt mit Hilfe all- 
gemeiner topologisch-funktionalanalytischer Betrachtungen in weitreichenden 
Fallen zu Existenzaussagen (iiber Existenz eines Fixpunktes), so zum Beispiel 
bei nichtlinearen Randwertaufgaben elliptischer Differentialgleichungen bei 
SCHAUDER und LErRay [20], [12], bei nichtlinearen Integralgleichungen bei 
RorueE [18] und andern. 

3. Bei linearen Aufgaben kommt es haufig vor, dass die Aufgabe zu einer 
Klasse von Aufgaben gehort, fiir die die Existenz einer Lésung durch besondere 
Sitze gesichert ist. (Das trifft auch fiir einige Klassen nichtlinearer Aufgaben zu.) 


9. Gleichungssysteme monotoner Art. Einschliessung 
bei der Relaxation 


In dem Gleichungssystem 
Ax=r (75) 


seien die gegebene nichtsingulare Matrix A = (a;,) und die gegebene Spalten- 
matrix 7 = (7;) reell; x = (x,) ist die Spaltenmatrix der Unbekannten; A ist 
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zugleich der Operator T. Die Bestimmung von «x ist eine Aufgabe von mono- 
toner Art, wenn ¢ = 0 aus Ae = 0 folgt, und das gilt genau dann, wenn alle 
Elemente von A~1 nichtnegativ sind. Dann werde A eine «Matrix von mono- 
toner Art» genannt. Ein einfaches, bei vielen Anwendungen erfiilltes, hin- 
reichendes Kriterium lautet: Die Matrix A ist von monotoner Art ((5], Arch. 
Math. 3, 373 [1952]), wenn gilt: 

1. Es ist a;; >0, a;, <0 fiir] + Rk. 

2a) Es ist das «schwache Zeilensummenkriterium » erfiillt 


n ZO frig ‘Ly, 2 0 | 
45 I 


(76) 
h=1 | > 0 fiir mindestens einen Wert von 7. | 


' 2b) Die Matrix zerfallt nicht, das heisst, es ist nicht mdéglich, die Indizes 
ees 7 SOM 07, 6-5 On; 07, -+-, On (mit 1 = m Sw — 1) zu numerieren, 
dass Q,0, = 0 fir y=1,...,m, w=1,...,%—m gilt. An Stelle von 2a) und 
2b) kann auch 

2c) das «gewohnliche Zeilensummenkriterium» treten: 


ORM COT ge Aes, (77) 
had 


Dieses Kriterium ist zum Beispiel erfiillt bei den Differenzengleichungen, die 
der ersten Randwertaufgabe der Potentialtheorie entsprechen, so dass man 
bei Anwendung der Relaxation bequem Schranken aufstellen kann (Zahlen- 
beispiel bei [5], Z. angew. Math. Mech. 32, 83 [1952]). Auch bei den Differen- 
zengleichungen bei anderen komplizierteren Randwertaufgaben und auch bei 
gewissen Typen nichtlinearer Gleichungssysteme liegt monotone Art und da- 
mit die Méglichkeit der Einschliessung vor [5]. 


10. Randwertaufgaben monotoner Art 


Es liegen wieder die Randwertaufgaben (51) und (52) vor; nun werde der 
Operator T durch 


Tv = Liv] — o(%;, ») (78) 
festgelegt, wobei der Definitionsbereich F aus den Funktionen v besteht, welche 
die Randbedingungen (52) erfiillen und partielle Ableitungen so hoher Ord- 


nung besitzen, wie es fiir U,[v] und L{v] gebraucht wird. Fiir zwei Funktionen 
v, w aus F gilt dann mite=v—w 


Tv—Tw=L{e]+eA(x,), (79) 
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wobei nach dem Taylorschen Satz 


Op 


Plt 2) — glx, w) = & (SE = A(x;) (80) 


ese 
(mit w als Zwischenstelle zwischen v und w) bei festem v und w die Funktion 
A(x;) eine Ortsfunktion ist. 

Jetzt sei H ein beziiglich wv konvexer Bereich im (x;, “)-Raum, der v, w 
und damit auch w enthalt. In H werde A als = 0 vorausgesetzt; man wahlt 
oft zweckmassig H so, dass H eine Lésung u von (51) und (52) und eine Nahe- 
rung v enthalt; mit w = u ist e dann der Fehler ¢ der Naherung v. 

Nun kann man in vielen Fallen auf Grund spezieller Eigenschaften der 
Aufgabe (51) und (52) zeigen, dass aus 


L[e]+ ¢A(x;) 2O0in B, Ufe]=Oaut I, (81) 


é = 0 bei nichtnegativen A(x;) folgt. Dann sind (51) und (52) von monotoner 
Art; es kommt also auf die Diskussion der linearen Aufgabe (81) mit homo- 
genen Randbedingungen an. So sind zum Beispiel folgende Klassen von Rand- 
wertaufgaben von monotoner Art: 

I. Die Aufgabe 


— [B(x) u(x)’ + F(x, w) = 0 (82) 
mit den Randbedingungen u(a) = u,, u(b) = m, oder 
w(a)— cua) =y_, u'(d) + dud) = yy (83) 


[p(x) na < x <b stetig differenzierbar und > 0, a,, Up, Ya, Y» Zegebene Kon- 
stanten, c= 0, d => 0] ist im Bereich H [a Sx Sb; u(x) Sus u,(x)] von 
monotoner Art, wenn dort 0//0u existiert und in H entweder 0//0u > 0 oder 
0f//Ou = 0 und dann c2 + d? > 0 ist. 

II. In (51) sei 


= 07 u = Ou : 
Lu] = ect crM Tr Dest wren aye (84) 
j,R=1 J ° q=1 J 
Dabei seien a;,, b;, c gegebene, in B + J” stetige Funktionen von Wa cn dy ee 


as 0 und die Matrix der a;, in B durchweg symmetrisch und positiv definit. 
L sei als Randflache von B eine abgeschlossene zusammenhingende (n—1)- 
dimensionale stiickweise glatte Hyperflache. Als Randbedingung sei entweder 


u=yaut I (85) 
oder 


Ou 
=~ —ou=daufl (86) 
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gegeben, wobei vy die innere Normale und y, @, 0 auf I” gegebene stetige Orts- 
funktionen mit 9 > 0 sind. Die Randwertaufgabe ist in einem beziiglich u 
konvexen Bereich H des (x,, u)-Raumes von monotoner Art, wenn dort Og/Ou 
existiert und S< 0 ist. 
Ith 
O4u o4 4 
Adu = 5 +255 x i 


=7(x, ¥) in B 
u=Y,(s), Au =y,(s) auf I, 


mit s als Bogenlange auf der Randkurve J" des ebenen, einfach zusammen- 
hangenden Bereichs B; y,, y,, y sind gegebene, etwa stetige Funktionen. Der 
Operator Tv = AAv—r(x, y) (dabei muss v wie oben die inhomogenen Rand- 
bedingungen erfiillen) ist dann von monotoner Art. 

Die Beweise fiir das Vorliegen der monotonen Art findet man in [5]; sie 
beruhen im wesentlichen auf dem Hilfssatz vom Randmaximum (vgl. zum 
Beispiel [5], Z. angew. Math. Mech. 32, 203 [1952]): 

Fiiv eine nichtkonstante, mit stetigen partiellen Ableitungen bis zur zweiten 
Ordnung einschliesslich versehene Funktion w(x1, ..., X») set bei dem Differential- 
ausdruck (84) unter den in II genannten Voraussetzungen Liw| <0 in B. Der 
Grosstwert von win B+ I" sei M. Dieses Maximum wird im Falle c = 0 stets 
und im Falle c= 0, wenn iiberdies M = 0 1st, nur auf dem Rande I ange- 
nommen. 


11. Rand- und Gebietsmethoden 


In Nr.10 wird stets gefordert, dass die Naherungsfunktion v die Randbe- 
dingungen (manchmal die homogenen, manchmal die inhomogenen) erfiillt. 
Die hier gegebene Abschatzungsart eignet sich daher besonders fiir die gewohn- 
lichen Differentialgleichungen, bei denen man in der Regel beim Naherungs- 
ansatz die Randbedingungen erfiillen wird, und bei den partiellen Differential- 
gleichungen fiir Ansitze, die die Randbedingungen befriedigen (sogenannte 
Gebietsmethoden, wie zum Beispiel beim Ritzschen Verfahren). 

Es sei zum Beispiel die Randwertaufgabe linear 


Ll eer he. hn) (87) 


wobei L[u] durch (84) gegeben ist und die Koeffizienten die dort genannten 
Voraussetzungen erfiillen, und die Randbedingung laute 

Ulu] =y auf I (88) 
und sei die Bedingung (85) oder (86) (ebenfalls mit den dort genannten Voraus- 


setzungen). 
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Aus 


Lp | a0; U[z|=90 
folgt dann z <0 und aus 
|Z[1]| S L[+2], U[4] = U[z.] = 0 


folet somit |z,| Sz, da man nur den letzten Schluss auf z= 2, — z, und 
Z = — 2, — %, anzuwenden braucht. 
Nun sei 2(x;) eine Funktion mit 


Uz) =O aul yr zea Soins (89) 


wobei A eine Konstante ist; eine solche Funktion wird sich in vielen Fallen 
leicht angeben lassen. Ist nun v eine (88) erfiillende Naherung, so gilt fiir den 
Fehler € =v — u: 


L[t]=LlvJ—7, U[¢]=9, 


also 


[5] = |v — «| < 4 Max |r — Zfe]]. (90) 


Bei partiellen Differentialgleichungen arbeitet man haufig auch mit An- 
sdtzen, die die Differentialgleichung, aber nicht die Randbedingungen erfiillen 
und bei denen man etwaige freie Parameter zur méglichst guten Anpassung 
an die Randbedingungen verwendet (sogenannte Randmethoden). In diesen 
Fallen kommt man oft mit dem am Schlusse von Nr.10 genannten Hilfssatz 
unmittelbar zum Ziel (vgl. [5]), was hier in diesem Zusammenhang genannt sei: 

Bei der Randwertaufgabe (87), wobei wieder L{w] durch (84) gegeben ist 
und die Koeffizienten die dort genannten Voraussetzungen erfiillen, und der 
Randbedingung (85) sei v eine die Differentialgleichung L{v] = 7 erfiillende 
DaDee heaton Von der Fehlerfunktion w= v— wu sind die Randwerte 


w= V—u=v—y bestimmbar. Nimmt w den Wert Null an, so gilt im ganzen 
Cablet B 


Wmin = We cag) (91) 


nimmt w nicht den Wert Null an, hat also w ein festes Vorzeichen, so hat auch 
win B+ I’ dasselbe Vorzeichen, und fiir den Betrag der Fehlerfunktion gilt: 


[27] S | |mae- (92) 
Im Falle c = 0 gilt (91) unabhangig davon, ob w den Wert Null annimmt oder 


nicht ; dabei bedeuten Win, @max Und | © | max das Minimum bzw. Maximum 
von w oder |w| auf dem Rande J", also bekannte Gréssen. 
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Die gleiche Aufgabe der Einschliessung der Lésung wu in Schranken behan- 
deln Diaz, GREENBERG, SYNGE und CoopERMAN [7], [8], [10], [23] mit einer 
funktionalanalytischen «Hypercircle method» und verwandten Methoden. 


Beispiele 


I. Nichtlineare Differentialgleichung: 


iL 
Au = u,, 4 pane ee Ta a ot 1 


u = 1am Rande x? + y? = 1 (stationdre Temperaturverteilung u(x, ) in einer 
Kreisplatte ; am Rande konstante Temperatur wv = 1, im Innern Warmeabgabe 
proportional ~ + u?). Es ist leicht, die Randbedingung erfiillende Funktionen 
Up, Uy mit Au, = uy) + uz anzugeben; es soll hier nur ein ganz grober Ansatz 
durchgefiihrt werden. Fiir 


“,=1+ (1—7?) (« — =+ 5% 7) wird Au =—4 (a+ = 


und uw, berechnet sich zu 
1 eet 
Uy = 1 sie: 4 Au). 


1. Existenz einer Lésung: Dazu geniigt es hier, die Uberschlagsformel (67) 
zu benutzen. Eine das Vorzeichen nicht wechselnde Eigenfunktion ist 


Wey) = Jo (Var) [Bessel-Funktion Jy; A= A, ~ (2,4048)?]. 


Der Bereich F sei die Gesamtheit der in x?+ y? <1 stetigen Funktionen 
u(x,y) mit Ow <1; dann ist Nigg=(1+2 4) maz=3 und nach (67) K=3/A; 
K/(1— K) = 1,07. Nun werde « zunachst etwa so gewahlt, dass u, und up fiir 
y = 0 iibereinstimmen; es wird u, = u, = 2/3 fiir « = — 1/3; man berechnet 
dann up, “, nach (63): 


|| >) — Uy || = Obere Grenze dienes 0,15, 


und nach (13) die «Kugel» S: 
uw — u| $0,161 oder |w— | <0,161- jy (V/47). 


Diese Kugel S gehért zu F, es ist also die Existenz einer Losung in F gesichert, 


die sogar in S liegt. : : : 
2. Untere Schranke der Lésung: Die durch S gegebene Fehlerabschatzung 
liefert zum Beispiel fiir den Mittelpunkt (0, 0) = (2/3) + 0,161 und ist noch 
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sehr grob. Man kénnte durch sorgfaltigere Wahl des Ansatzes die Fehler- 
schranken des Iterationsverfahrens verbessern; aber da hier die Aufgabe nach 
Nr.10 von monotoner Art ist, ist es viel bequemer, zur Abschatzung das Ein- 
schliessungsprinzip (73) und (74) zu benutzen und das Iterationsverfahren nur 
zum Existenznachweis. Da fiir den Wert « = — 1/3 bereits u, < up ausfallt, 
gilt beim Operator Tv = —dv + v + v? auch 


Tu, = —Mu, + 4,+ 4 S—-Au,+u,+u=0, 


also wegen Tu, <0 nach (74) wu, Su. 

3. Obere Schranke der Lésung: Es wird dazu nur noch eine Funktion u, 
mit u, = u, bendtigt. 

Fiir u, = 1 — (1 — 7?) [0,32 + 0,05 r? + 0,006 74] wird u, = uw) und damit 
u, =u; mithin ist « in Schranken eingeschlossen; speziell fiir y= 0 wird 
0,6666 < u(0, 0) < 0,6800, fiir 7 = 1/2 wird 0,7396 < u(1/2, 0) S 0,7503. Natiir- 
lich ist es leicht, durch genauere Ansatze engere Schranken zu erhalten. 

II. Eine vielbehandelte lineare Aufgabe: Torsionsproblem fiir ein Quadrat 
B als Querschnitt. B sei in einer (x, y)-Ebene der Bereich |x| <1, |y| <1, 
und J’ der Rand des Quadrats; dann ist gegeben Au = 0 in B, u = (x? + y?)/4 
auf J’; es mégen harmonische Naherungsfunktionen v verwendet werden 


p 
“sv o, av, mit v,=Re(x+17y)*. 
= 


Die a, bestimmt man so, dass sich v den Randwerten # = (x? + y?)/4 anpasst, 
was sich zum Beispiel graphisch sehr rasch durchfiihren lasst. Benutzt man 
nur vy und 2, so unterscheidet sich (1/160) (47 — 8 v,) auf J’ von w um héch- 
stens 1/160, also gilt nach (91) in ganz B 


VFELS oy) SU 


1 
“— Te ( 


zum Beispiel im Mittelpunkt 0,2875 < u(0, 0) < 0,3000. Nimmt man noch Up 
hinzu, so kann man 


1 
v= > (1,1786 — 0,1801 v, + 0,006 2.) 


verwenden ; auf /"gilt dann |v— «|< 0,00125, daher in ganz B |v—u|< 0,00125, 
speziell im Mittelpunkt 0,2934 < u(0, 0) S 0,2959. 


12. Das Newtonsche Verfahren bei nichtlinearen Randwertaufgaben 


Das Newtonsche Verfahren liefert bei manchen Typen von Randwertauf- 
gaben die Méglichkeit, zugleich die Existenz von Lésungen zu sichern und die 
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Lésungen in Schranken einzuschliessen. Es liegen wieder die Aufgaben (51) und 
(52) vor. Ausgehend von einer Funktion u(x,;) bestimmt man u, = % + 1% 
aus der linearisierten Randwertaufgabe: 


L[uq] = Llu] + Lino] = (x5. Mo) + No PulXs, Mo) in B, 


Ulwl=y, auf Il, (w=1,...,%), | (93) 


wobei ,, die Ableitung von m nach w bezeichnet. Allgemein ermittelt man 
beim gewoéhnlichen Newtonschen Verfahren die Folge der Iterierten u,,(x;) mit 
Uni = Un + Nn ach 


Linn] — Nn PulXj, Un) | P(%5, Un) Te L[u,,] in B, U, [Mn+] Sen: aut ie , (94) 
und beim vereinfachten Newtonschen Verfahren nach 
Lal — Nn Pul%s, Mo) = V(%4, Un) — L[uplin B, Ulenal=y, auf L,. (95) 


Die beiden Arten unterscheiden sich also nur in dem Argument bei 9, (vel. 
Nr. 4). Die Transformation g = T/ wird fiir eine stetige Funktion /(x;) durch 


Le] = (xf) + @—f) eults %) in B, Ulgl=y,auil, (96) 


festgelegt, wobei wieder vorausgesetzt werde, dass die lineare Randwertauf- 
gabe fiir eine Funktion w 


Llw] — w pul%j, %) = 7, U fw] =Oaut I, (97) 


fiir beliebiges stetiges 7 eindeutig, eventuell mit Hilfe einer Greenschen Funk- 
tion G(x,, €;) nach (54) lésbar sei. Fiir zwei Funktionen g; = Tf; (j = 1, 2) gilt 
dann mit Benutzung des Taylorschen Satzes: 


Le, — Gel — (21 — Sa) Puls, Mo) = P(%i fh) — (%H fe) — (fr — fe) Pul¥ i Mo) 
= (fy — te) [pul f) ~ Gals 0] 
= (fa — fa) f= 40) Purl te); 
O81 = 20] =—0. 
Dabei ist } eine Zwischenstelle zwischen /, und /, und Uy eine Zwischenstelle 


zwischen f und uw). Dann wird 


q 
i 


{git Sp = | G(x, €4) [f1(&;) a fo(E5) | [/(&:) = u9(&;).| Puul $i, U9 (&;)} dé; . 
B 
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Mit der Norm (63) folgt (64) mit der Lipschitz-Konstante 


K = Obere Grenze way | |G(x3, &5) | W(Es) Puulés, Mo(;)] 48; | 
in B (%;) 
B 


x Max | f(x) — Uo(x;)| - | 
in B 


Grundsatzlich ist A < 1 erreichbar, sofern man erstens “, hinreichend nahe 
bei einer Lésung gewahlt hat, das heisst, sofern u, — u) und damit auch f — up 
geniigend klein ausfallt, und zweitens eine Greensche Funktion existiert, das 
heisst, die zu (97) fiir y = 0 gehorige Eigenwertaufgabe keine Eigenlésungen 
besitzt, ein Ausnahmefall, der durch kleine Abanderung von w#, vermieden 
werden kann. Uber Méglichkeiten, die Greensche Funktion bei der numerischen 
Rechnung zu umgehen, und ein Zahlenbeispiel vergleiche [6]. 


13. Anfangswertaufgaben 
bei hyperbolischen Differentialgleichungssystemen 


MORGENSTERN [15] behandelt das System von Differentialgleichungen fiir 
zwei Funktionen u(x, y), v(x, y): 


Ou ov 

jy Aut bv+e, a9. ne te eae 
wobei a, ..., f gegebene Funktionen von x, y sind und auf der Anfangskurve 
x + y = 0 die Werte Null fiir w und v vorgegeben sind. Es wird ein Operator T 
definiert, der einem Funktionenpaar {u, v} das Paar T{u, v} zuordnet: 


T{u, v} 


V 


| ‘i (au+tbu+c) dé, / (du+ev-+ f) dn. 


H=—-x 


—_—_ 


Als Norm eines Funktionenpaares kann man 


| {, v}|] = Max (en |o(x, 1) 


inB \W(x,y)’ W(x, 9) 


verwenden, wobei B etwa das Dreieck x < ByS-a,x*x+y2>0mita< B 


sei und W(x, y) in B eine festgewahlte stetige positive Funktion ist. Fir zwei 
Funktionenpaare w,, v, und up, v2 gilt dann 


x 
/ a(uU, — Uy) + b(v, — v4) dé 


|| Ti, v1} - T{tig, vq}. = Max a Wiss wees 


SK |{u, — uy, 0, — Us} | 
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merit 


K S Max 


Ma ([a(g, 9) | + [OE »)) WE, ») a 


, 


1 
W(x, 9) 


: | 

| 
€=-y | 
wircse | fae es nb. nan | 


MORGENSTERN setzt W(x, y) = e*(**”), dann folgt 


—< | 5 
K = Max ( T - ae iy 


wahlt man also L > Max (|| + |d|, |2| + |e]), so wird K <1. 


14. Integralgleichungen 


Bei der linearen Fredholmschen Integralgleichung zweiter Art 


u(s) —A | G(s, t) u(t) dt = h(s) 
e.g a 


mit gegebenen (etwa stetigen) A(s), G(s, 4), gegebenem A und B als Integra- 
tionsbereich bei einer oder mehreren unabhangigen Veranderlichen sei der 


Operator g = Tf durch 
+4 | Gs, 1) HO at 
B 
festgelegt. (Allgemeinere Operatoren betrachtet SCHMEIDLER [21].) Dann gilt 
bei der Norm (63) fiir zwei Funktionen g; = Tf; (j = 1, 2) wieder (64) mit 
= {A| Obere Grenze 4p oy { I (s, )| W(t) 


BUCKNER ((3], [4], S. 68ff.) gibt eine Verallgemeinerung des Iterationsver- 
fahrens 


tensa(8) = Puan Ha(s) + ACL — Pos Baas) | Gl G(s, #) tg(l) dt + (1 — Ings) MS) 


an, bei welcher die Folge der ,(s), sofern 2 kein Eigenwert ist, bei passender 
Wahl der #,, stets konvergiert. 
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KANTOROVITCH ([11], S. 90ff.) untersucht auch Systeme von Integralglei- 
chungen und gewisse Typen nichtlinearer Integralgleichungen; LONSETH [13] 
betrachtet Integralgleichungen erster Art, PICONE und FIcHERA [17], [9] neben 
diesen auch Integrodifferentialgleichungen. 
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Summary 


In this collection of results a formulation of the basic theorem of iteration 
methods is given, which is particularly suitable for practical applications since it 
permits testing the presumptions in a straight forward way. A uniform existence 
and error estimation theory results, which is valid for many various types of 
ordinary and partial differential equations. The well-known criteria for iteration 
methods for linear equations are automatically included and are extended to 
apply to systems of non-linear equations. Initial and boundary value problems 
for ordinary and partial differential equations are also treated. 

For numerical calculations and especially for non-linear equations, it is rather 
important to chose a suitable norm. In problems of a monotone nature, a check 
on the accuracy may be achieved by the simultaneous application of the theorems 
on iteration methods and on error estimation. 


(Eingegangen: 11. Juni 1953.) 
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1 
Transonic Flow Past a Wedge at Zero Angle of Attack’) 
By Leon TrILt1NnG, Cambridge, Mass., U.S.A.*) 
While recent theoretical investigations of transonic flow past wedges, by 


J. D. Cote [1]’), GuDERLEY and YOSHIHARA[2], MACKIE and Pack[3], VIN- 
CENTI and VAN WAGONER(4] and experiments by LiEpMANN and Bryson [5], [6] 


Fig. 1 
Transonic flow past a symmetric wedge in physical plane. 


and GRIFFITH [7], {8} have gone far to explain the structure of certain transonic 
flow fields, there is no published discussion of the pattern which results from 
a local supersonic region embedded in a subsonic field. The simplest flow of that 
type is past a thin symmetric closed wedge at zero angle of attack. 

Let the chord of such a wedge by 2c, its thickness 2 ¢, and its semiwedge 
angle 6 so that 6 = t/e-<1. The free stream Mach number M is just below 


1) This investigation was carried out under contract AF33-(038)-22184 with the United States 
Air Force. 

*) Massachusetts Institute of Technology, Cambridge, U.S.A. 

3) Numbers in brackets refer to References, page 375. 
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unity so that 1 - M? <1. The sonic line and a shock JG (Figures 1 and 2) then 
bound a finite supersonic region CIG as can be seen on interferograms [5], [6], [12]. 
Disturbances downstream of this region affect the supersonic flow through their 
influence on the upstream subsonic flow, so that the influence of the subsonic 
flow cannot be limited to part of the supersonic field. 

If the velocity perturbations w, v from sonic velocity a* are small (u, v/a* <1), 
the strength of the shocks DE, FG is of order (wu, v/a*), and the entropy increase 


Fig. 2 
Hodograph of transonic flow past a symmetric wedge. 


across them is of order (w, v/a*)*; the flow is therefore assumed to be isentropic 
and irrotational, and the transonic perturbation theory may be applied. In 
terms of the dimensionless velocity perturbations: 

u(y +1) pba set 


uo = as ; ar 
the equations of motion are: 
, Ou’ ov’ 0 (1) 
eroy Ov ; 
Ou iy 0. (2) 
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The Jacobian of the hodograph transformation 1s 


O(u’, v’) it 08 Neate ai 3 
ere ict 3) 
so that the transformation is regular and order-reversing in the subsonic domain, 
but may have singular lines in the supersonic domain. The equations of motion 
in the hodograph plane are: 


Oy Ox 
Pim te eee a) 4 
e Ov’ Ou’ ( ) 

Ov Ox 
pe A ae te 5 
Ou’ ra hig ( ) 


Equations (4) and (5) are elliptic in the subsonic region (w’ < 0) and hyperbolic 
in the supersonic region (w#’ > 0), where their characteristics are defined by: 


hae 
v=ctsw™. (6) 


The following boundary conditions are satisfied by the functions x, y(w’,v’), 
with the origin in the physical plane at the center of the wedge. Far from the 
wedge, the flow is subsonic and undisturbed, so that if w; denotes the free 


stream velocity, 


: Cc 

lim x, y(#, v) = ————__— (7) 
ul Vuy— 
v0 


as in an incompressible fluid. 

This boundary condition implies that the velocity is identical far upstream 
and far downstream of the wedge. No inconsistencies arise since the drag is of 
second order, and second order velocity perturbations are neglected. 

Along the streamline y = 0, from A to B (Figures 1 and 2), the velocity is 
parallel to the free stream and decreases from the free-stream value u; to the 
stagnation value u’ + —oo. From B to C, the v’ velocity takes the value 


v =u =(y+1)6 


while the w’ velocity increases from —co at B to zero at the corner C. At C, 
the flow turns the corner CD through a Prandtl-Meyer expansion CD. 

Because the velocity at the end of a Prandtl-Meyer expansion depends only 
on the angle of turn, while along EF it depends on the free stream Mach number, 
a weak oblique shock may be necessary to adjust conditions at the corner: the 
flow overexpands within the supersonic region and is recompressed through an 
oblique shock, wich is sufficiently well represented by an isentrope. The jump 
takes place from one sheet to another in the hodograph (Figure 2) similar to 
those made in supersonic wing theory. 
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From D, the supersonic flow decelerates to F, where a locally normal shock 
FG closes the supersonic region (9). From G, the subsonic flow decelerates to a 
rear stagnation point H where it turns back to free stream direction; it finally 
accelerates to free stream velocity wu, and the line y = 0 is closed at A. The 
variable x may reach the value c before P. The flow must then turn back to 
v’= 0 through a weak oblique shock followed by a curved shock behind the 
airfoil. This more complicated case is not discussed here. 

The function x(u’, v’) goes from — co to —c along the upper edge of the 
cut AB; between B and C, it increases from —c to zero; it remains constant 
along the characteristic CDE, then x increases from EF to F, remains constant 
across the shock FG, increases until it reaches the value +c at the rear stagna- 
tion point H, and increases from ¢ to + co as it moves from H back to A. When 
y(u', v') is known, x(u’, v’) is determined from (4), (5) up to an arbitrary 
constant which adjusts the position of the wedge along the x-axis. The bound- 
ary conditions across the shock FG are integral relations 

G G 


[vau= | xdu=0. 


F F 


The condition y = 0 between F and G is the simplest particular solution of this 
integral relation. 
To sum up, the following boundary conditions are satisfied by x, y(u, v): 


x,y have a specified singularity at A. 


x takes the values --c as u’ > —oo, v ¥ 0. 
On the subsonic contour CBBAHHF’, y vanishes. 


On the characteristic CC’D, x and y vanish. 


| 8) 
| ( 


The solution in the supersonic region F’ C’ DEF’ may be found separately by 
a characteristic calculation with boundary conditions specified along C’F” by 
the previous solution and with x, y vanishing on C’D, and y vanishing on EF’. 


Construction of the Solution in Region CBBAHHF’C’ 


The solution is determined in several steps. First, the appropriate singular 
solution and stagnation solution are constructed. Boundary conditions on BC ; 
F'G are satisfied by a system of images. Then, the solution is made to vanish 
approximately along CD. Finally, it is adjusted along the cut BAH. 
It is convenient to introduce the following new variables 


wee ; * phys i A pe ; 
yy =e C V phys»? Mie 6 2 U6 f . 2 uy’ 3/2 
es 
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The lines v’ = + v) then become v = + 1 and the flow is entirely described by 
the parameter 


2) MELO een 5 
eee = a 
ay 3 (y 43 1) } 3 ( 7) , ( ) 
where &j is the transonic similarity parameter. 

By separation of variables, one obtains the component solutions 


ys GEA 2S CyglA2) 9 ot user att Gg eels (10) 


where 4 is any real or complex number and C, represents any linear combina- 
tion of Bessel functions of order -+ ». | 
Another type of solution is found by introducing the variables 


te w+ By)? (11) 
and seeking solutions of the form 
Vara" f(t), ee See gee (12) 


Tomotika and TAMADA[10] point out that ¢ can be made equal to unity at 
any point of the z- (or negative w-) axis by the choice of an appropriate imagi- — 
nary value of B. In particular, if B =72,, ¢ = 1 at the point (0, w,) of the hodo- 
graph plane which corresponds to the point at infinity in the physical plane. 
The functions /,,(t), g(t) satisfy the equations: 


t(1—?%) fh -[Fm-y+ oq n+" of Deatly 


(13) 


so that two linearly independent solutions with singularities at ¢ = 1 are given by: 


yy ¢ 
y, =u" (—2)n"8 F{ Zn +S. nog t), 


6 ee 
l4a 
MY = yo” (= fern F{ 2% -bew Wit 1 4 “7 ( ) 
2 6 , 3 zr), 
a. ? 2m-+ 3 m 1 
Cs m (__4\m/3 — Fiz) 
EE te eT ARIE, 
(14b) 


= Ie 2m 7 m 2 5 
ry =u (44218 ——, _ ; 534) 
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These are related to the solutions (10) by the relations: 


ie.) 


oie z Ree eb aN 
BP fe Morin) F(A) AY dd = 28 (v4 i nytt hap 
0 
34-1 3yu+2 2 — 2 
Dares ;—: 
6” Sa eke Zz)? )? 
( (ui + 0 2;) (15) 
Pair : a TU3 w+ 1)/3) 
A(vu+t2) t—1 <2 1/3 (1+ 3 w)/3 
ae |e \ Digan {dh = 2a (gy 2) Us 8a) coe. 
0 
3u4+1 Ose Ge ah — 2 
x F Grima 6 nosey 
ser 8 e I[(3 w — 2)/3] 
—A(v+izy = 2/3 3 u)/3 
BP Raga ah 2 (+ iaye-sme 
0 
3u— 2 Snes ih ah — 2% 
wal @ ? 6 ae Tas, feat (16) 
fe ee hay A aA 2 tg) P(2 Glee 


0 
PES jp 3 ps 


5 a ge 
6 4 6 ap ae are 


x F( 


Substitution of the integrals (15), (16) into (4) and (5) shows that a solution 
of the system is obtained by taking 


eT A 


When the Fourier-Bessel integrals exist, they are convenient to satisfy 
analytic boundary conditions along lines v = const; the hypergeometric func- 
tions are useful in finding the singular parts of the solution. 

It is easily verified that 


lim (1 — #) = lim > (w— 1), (17) 


Uy 
t—1 U—> Uy 


so that the singular part of the solution behaves as (l — t)~ 1? near bm ithe 
singular part of the solution y,(v, v) is an even function of »% and x, is an odd 
function since the flow far from the wedge is similar to an incompressible flow 
y, is therefore proportional to the real part of (1 — ¢) 5 and x, to its imaginary 
part. It follows from (17) that the proper singularity is found by selecting 
n=1, m=1/2. After separation of real and imaginary parts, the singular 
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solution becomes 


fo,e) 


ee a a 


0 


x fara AL [Tig(A 2) + J_1jg(A 2)] — cosa 21 [Jzj(A 2) — J_1/s(A ‘| day (8a) 


=Re{—(3)" e Fa F(-2,555 ;1-4I, 


lo, <) 
EIR yh 
ce ed ere 
0 


< [cosa 24 [J_9j3(A 2) + Jojg(A2)] = aes [J-2j5(4-2) — Jojg(A ‘| da ; (18b) 


= m{ (3) east mee 


When v = 0, the integrals (18a) and (18b) are evaluated as follows (see [11], 
page 405) : 


pane z 1/3 z —1/3 
Tae eoAlerar cs ah sr Cag 
Vs =U, (z > 2) 


(19) 
%,=0, (2 <4) 


%,= + (5)" = - jeos (3 Sin 4) oa yn (= sin-} 2)]. (z > 2,) 


This solution exhibits the required singularity and symmetry properties, but 
does not satisfy the requirements of a stagnation point, since lim x, = 0(z~ 8) 
instead of a constant value x,, = -+ c. To find the stagnation function, consider 
the solution 


Bees 4/9) ogee es Nata Re 
A eae I'(5/3) eres: | FS. 3) 3) ware (20) 


As |2| 00, X,,—> (2-4? I'(2/3)] since the hypergeometric function is finite 
at unity. This result holds for all values of v, as it should since at the stagnation 
point, v takes all values in the range 0 <|v| <1. In order to obtain the ap- 
propriate antisymmetry across the cut, one defines: 


2) . 1 EB) 
1 — ettria 2 4 VE 1 
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Then, %,;, the imaginary part of X,,, has the proper behavior at infinity. 
On the other hand, when z < 2,, x,, vanishes. In terms of Fourier-Bessel 
integrals: 


= ies 2/% een il ee 5 
ais Fay ? ie anes [e CsA 27 = Ve sind a| Joe) dA (21) 
0 


[o,e) 


) pew C cos A.2, = ee sind a| Faajg(4 2) Ada. (22) 


4\5/6 21/8 
T'(2/3 
0 


st (+ 


In particular, along the z-axis: 


= ni ees) Ala 
a / 
ov. Re "8 Beas : = : (7 aye] (23) 
When z > 2, one obtains 0x/0v = 0y/Ou = 0, so that 
oy hy EU Oe ee 
Ou VT'(2/3) ~ iz “Ga gyi (2 > 2;) | 
24 
merle (7/6) “Als eo (24) 
— (2/3) U 1G . (22 as Z3)716 * (z 21) 


The symmetry conditions are satisfied, but y,, does not vanish along the cut. 
To cancel the contribution of y, + y,, on v = + 1, one now puts two sin- 
gularities of the same strength and opposite sign at v= + 2+ 72,; to cancel 
the effect of these first correction solutions one requires solutions of identical 
strength and sign as y, + y,,at v = +4-+74,, and so on. 
In terms of Bessel-Fourier integrals, the image correction is written as 
follows: 


(Vs =a Vst)e =| Did, z) je se (ae en7thn (ent ee “9 dA 
6 1 


cosh 4 


= [ 00, 2) le” e¢ con] dh (25a) 
0 


sinh A (1 — |v]) 
= | 00, 2) cosh A = di, 
0 
where (A, z) represents the combination of Bessel functions under the integral 
* sign in (18), (20). Conditions along v = + 1 are satisfied. Conditions along OA 
are also satisfied because of the symmetry of the images. However, y does not 
vanish along the cut BAH. 
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The corresponding relation for x(w, v) is 


eae _, sinhdv 
(x, + ta)o= | PU z) [+ en 7ll — e-4 | dh 


(25) 


~ + [w, 3) OS ea 


where Y/(A, z) corresponds to ®(A, z), and the region of convergence is the same 
as for y. 
One is then led to investigate solutions of the form 


le 2) 
ya Aye + Va. + PTE AES z)[a, cosnav+ bd, sinnz | 


cosh iv ~ sinh A (1 — |v|) 26 
" anf cosh A A) + suid F(a)| 78) 


x [Tijs(A z) + Ja1)3(4 z)| dd , 


1/3 
% = A Xe + Xou, + 28 (5 ) SI Kyy(n 2) [b, cosmav+ a, sinn av] 


ame crap Ape G(A) cosh A ( (1 — |v}) F(A)| (27) 


cosh A sinh A 


X [J_29/3(4 2) — Jojg(A 2)] da 


where (x, y),, represents the singular solution, (x, y),;, the stagnation solution 
both corrected to satisfy conditions on v = ++ 1. The series is needed to satisfy 
the condition on CC’. The function G(A) is used to adjust conditions along the 
lines v = + 1 and the function F(A) to adjust conditions along the z-axis when 
all other corrections have been made. 

To satisfy the condition y= 0 along CC’, following GUDERLEY [2], one 
replaces it by the equivalent conditions 


y=0onCN, y=~% = 2% _ on NC’ 


since these conditions insure the vanishing of the solution in the entire region * 
CNC’ and in particular along CC’. In the analysis, the Bessel functions J 41/3 
are replaced by their pe gs) SA values. 
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Along NC’, one therefore has to satisfy the conditions: 


Dy %m CosN 7+ By sinn av = flv) , 
: (28) 
De; (x, cosnau+ B, sinn 2 v) = g(v) . | 


Here 
An = An [Jijs( mt) + J 1)3(” ze) |? ay Vn , 
Pa = oF AAG It) te J_1)3(” ze) Be by, Vn , (28a) 


4 a, \aidu) [Cr Just 6) F Ot" J aalm aw 6) 
rm 67 [Sara(m mo) + J_ara(n x €)) 


and 
/(2) ae =A Vs, > Vete)em1 


3 ff (CO. G sinh A (1 — v) f 
fi | ( coshi °* sinha F) K,)3(A) 4A , 


(28b) 


Sl Spree hay 


(s2884 v sinh A (1 — v) 
+r 2/( cosh A ae sinh A ) K39(4) dh . | 
Since /(v), g(v) are defined only in the half-range 0 < v < 1, it is convenient 
to continue /(v) as an odd function of v so that 


PIs (29) 

a 

1 
Ya = 2 | ft) sinn mv dr . (29a) 

0 
Writing now: 
g'(v) = g(v) — D/ An Ynsinn ay, (30a) 
1 


and differentiating (28) one obtains the equivalent system: 


1 


Be <1) 


ys An (On cos v + d, Sinn 7 v) = g’(v). 


1 


Dd, # (bn COSHH — Hy Sinn tv) = 0, | 
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When the function £(U) defined by: 


0 0 
Oe = [ u) Gosney Ud eee = | 60) snmxUdU (32) 
=a -1 


is substituted into (31) and the range of the integral and g’(v) are made the 
same by the transformation y = v — 1, the system (31) can be written 


0 
[0O) Seno eorey WU = 20"). (-1<y<0) (33) 
—1 


cosx U — cosa p 


This integral has meaning if the Cauchy principal value is taken when y = U. 
To solve this integral equation, the angles U, y are replaced by angles which go 
through a full period in the range of integration: 


y=%, v=%, gty)=%@), oU)=7w), a=. 


= 


It can then be shown by expanding # and g” into Fourier series and carrying 
out the integration, that: 


p! 


NM Oy = -| g’(v) cosnaudv (meven), nma,=0 (n odd), (34) 
0 


1 


nO = at g(v)sinnavudv (neven), 6, =0 (mn odd). (35) 
0 


Returning to the original coefficients a,, b,, one obtains, using the orthogonality 
properties once more : 


1 


a, = — ia / g(v) cos 2 v dv , ( even) | 
“0 
, (36) 
ee fe | [g(v) + » f(v)| sinn 2 v dr, (n even) 
0 


where f, g are given by equation (28). 
Along v = + 1, y vanishes if G(A) satisfies the following integral equation: 


oO 


[GARD 4+ Fig] dd= Say Ky Qnn2) 87) 


0 
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or writing K,)3(2 1 x z) as a Fourier-Bessel integral ([11], page 424): 


[o.e) 


A (ae i Sy Aen Ea [Tra(A 2h Jo 1jg(4 2)|dA=0. (38) 


0 


This result must hold for all values of z (real and complex); therefore the 
integrand must vanish: 


Azn 4 / 
=-» A? + 4 n? 7? * (39) 


Along v = 0, x is continuous along OA and y vanishes along AB, AH when F(A) 
satisfies the integral equation set: 


J Acoth a F(a) [Jys(42) + Jg(42)] da =0, (<x) 
[ie (A) — (1 — sech A) G(A)] [Jiyg(A 2) + J ajg(4 2)] aA (40) 
= = yoi— fe? BA, 2) da. (2> 2x) 


0 


This diffraction type integral equation is not easily solved. However, since 
F(A) represents a correction term, equations (40) are approximated by 


Asin dA 2, (41) 


Tt 
OL SS era Tay 


When (39) and (41) are substituted into (36), one obtains an infinite set of 
simultaneous equations for the coefficients a,. But, since the a, decrease ra- 
pidly and the integrand decreases exponentially, it is sufficiently accurate to 
replace (39) by P 
a 
G(A) = — sepa: (42) 


The integrals in (36) are then estimated as follows: 


i 2 a 'coshA n a — 1 an Anasind anes 
n~ Vn sinha n x 2° Lt Ana)s? 
= (43) 
A as tanhéna\ A* Kor3(4 ” 2) dj. | 
eo fe 4/n2 : nN ie 1 + A? 


As na 1, because of the rapid convergence of the K); function, the asympto- 


ZAMP IV/23 
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tic value of Ky)3(A m2) as Ana oo are used and of the rest of the integrand 
asA>0O: 


re 3/2 rs ” oe Site —0-0511 
I, & 5 ae Ay nif Jie er FAO GA ox a (43b) 
0 
Similarly: 
a fy : mA [A Keys(4) — 0% Kyg(A)] 

In= = / (FQ) + G(A) (sech 4 — 1)] =4#"““—8 os ia | 
0 + (44) 

== 0.28608, | 


tLe — 
Vn 


The formulas (43) and (44) give the correction terms to be added in order to 


obtain the complete coefficients a,, b,. The constant A is adjusted so that 
x(0, 1) vanishes. The solution in the supersonic region may then be determined 
by the method of characteristics. 


Analysis of Results 


By means of the solution constructed above, several significant features of 
the flow are computed. They include the size and shape of the supersonic region 
CIG (Figure 3) and its variation with transonic similarity parameter (Figure 4), 
the variation of normal shock strength at the wedge surface with similarity 
parameter (Figure 5) the actual pressure distribution (a typical distribution is 
shown on Figure 6) and the variation of drag coefficient with similarity para- 
meter (Figure 7). 

Since y grows monotonically allong CD, in the subsonic domain (9), the 
position of point J, the apex of the supersonic region, is given by the condition: 
dy(0, v)/dv = 0. This occurs at the point: 


. co 
ty = (Z)" ——__ 7 ee, (45) 
Chins eI) pi eos 
21/8 4 288 (4/3) 28 ee 
Yo= aya TF 3878 T(5/3) 238 + 3172 Q18 T'(2/3) pa mus (455) 


As the transonic similarity parameter vanishes, yy grows beyond bounds. 
Since A remains bounded when z, approaches zero, the height y /c of the super- 
sonic region becomes inversely proportional to the similarity parameter. 

The width of the supersonic region is the distance from the corner to the 
intersection of the normal shock FG with the wedge surface. It is seen from a 
typical chordwise velocity distribution plot (Figure 6) that while x(w, +1) 
increases monotonically from =1 to zero as u increases from — oo to zero; 


— 
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Fig. 3.— Shape of supersonic region transonic flow past symmetric wedge. 
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Fig. 4.—Variation of size of supersonic region in transonic flow past symmetric wedge. 
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Fig. 5 


Strength of normal shock which ends supersonic zone at surface of wedge in transonic flow. 
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Fig. 6 
Velocity distribution at the surface of a wedge in transonic flow. 
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Fig. 7 


Generalized drag coefficient for symmetric wedge in transonic flow. 
© Total drag, x Drag of rear half of wedge, A Drag of front half of wedge. 


a(u, —1), increases to a maximum as & decreases from up to u = 0): 


(2 x) 18 atte aii pee 
© = “35 T(1/3) ier 96) 


As u continues to decrease, x decreases to a minimum G"; it then increases 
monotonically as w goes to negative infinity. In the region Xqn << % < Xp the 
function u(x, 0) is triple-valued. Therefore, a shock inust be placed at that 
value of x in the range xg” <x < xg at which the shock condition (xp) = —u(X¢) 

-is satisfied. The location and strength of the shock are automatically determined 
by the condition that w(x, 0) and more generally u(x, vy) be.single-valued. If the 
triple-valued region starts at x/c > 1, the flow must turn at H through a weak 
shock (or Mach wave) and become supersonic and parallel to the free stream 
behind the wedge. It becomes subsonic again by going through a curved shock 
symmetric with respect to the x-axis and concave in the downstream direction. 
The shock comes off the airfoil (x/c = 1 at F) for a transonic parameter value 
= £,= 0-12. For the range || < 0-12 the flow configuration is considerably 
more complicated than that investigated here. The position and strength of the 
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shock for —&, > 0-12 are plotted as functions of the transonic similarity para- 
meter &) on Figures 4 and 5. 

The pressure drag coefficient (Figure 7) is the integral of the pressure coeffi- 
cient component in the horizontal direction. Within the limits of the linearized 
theory, the pressure coefficient is: 


2(u — %) Lae ae 
ope = ape (47a) 


so that the drag coefficient is: 


0 1 
Gp = aaa (5)" am dx whe | | 
J 


od (47b) 


<< ete 


ON LIS Bee aes ie Si OY 
< | BD Gy" Eten + fleas fom 3 a, 
1 6 ig 

This drag may be considered made up of two related parts. The infinite 
series represents the drag due to the fact that the subsonic part of the flow-field 
is not symmetric force and aft. The integrals represent the suction on the rear 
of the airfoil, due to the supersonic region there. 

The drag may also be split into drag of the front of the wedge and drag of 
the rear. The drag of the front of the wedge agrees quite well with the drag 
computed by CoLE[1] and measured by Bryson{6] for the half-wedge. The 
drag of the rear is much larger, and due mainly to the suction on the rear of 
the airfoil. 


Conclusion 


The preceding pages describe a theoretical investigation of the steady plane 
flow of an inviscid ideal gas past a thin symmetric diamond wedge at zero 
angle of attack at speeds just below the speed of sound. A solution of the tran- 
sonic equation which approximately satisfies all the boundary conditions, is 
built up in the hodograph plane from an appropriate singular solution, a stagna- 
tion solution and a series of regular solutions. 

From this solution, a rather complete picture of the flow field is determined, 
including the variation of size and shape of the supersonic region, and the po- 
sition and strength of the shock wave which bounds it on the downstream side. 

The pressure drag is in good agreement with CoLr’s theory and Bryson’s 
measurements for the front of the wedge, and fairs in smoothly with VINCENTI’s 
curve at €) = (M = 1). The few available experiments on closed wedges give 
drags of the same order as the computed drags. 
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Résumé 


L’écoulement stationnaire plan d’un gaz idéal, sans frottement ni conduction 
de chaleur autour d’un profil mince en losange a l’incidence nulle, a été examiné 
théoriquement a des vitesses légérement inférieures a celle du son. 

Le probléme fut étudié dans le plan de l’hodographe ou la coordonnée (wu, v) 
est une solution de l’équation de TRICOMI, et ot x(u, v) satisfait 4 une équation 
semblable. Les conditions aux limites exigent que y(u, v) présente une singularité 
du type vy ~ (u, — u)-4/2, ot u, est la vitesse de 1’écoulement infini amont; y 
s’annule le long des droites v = + v9, qui représentent la surface du profil, et le 
long de la caractéristique descendante passant par (0, V9); les points d2 stagnation 
amont et aval sont définis par la condition lim #(—oo + U9) = F ©. 

Une solution a été élaborée qui satisfait approximativement a ces conditions 
et qui dépend uniquement, du paramétre de similitude transonique So: On en 
déduit que l’écoulement est subsonique jusqu’a l’aréte au milieu du profil, sonique 
sur l’aréte elle-méme, ou s’amorce une détente locale PRANDTL-MEYER limitée 
par un choc oblique faible. Suit une zone d’écoulement supersonique décéléré 
aboutissant a un choc normal situé soit sur la moitié arri¢re du profil soit dans 
l’écoulement aval. La résistance due a la pression, la position et lintensité des 
chocs, ainsi que la forme de la région supersonique, ont été calculées pour diffé- 


rentes valeurs de &. 


(Received: March 10, 1953.) 
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Verbesserung des Differenzenverfahrens von H. Gortler 
zut Berechnung laminarer Grenzschichten 


Von HERMANN WITTING, Freiburg i. Br.?) 


§ 1. Einleitende Ubersicht 


In der vorliegenden Arbeit wird ein Differenzenverfahren zur naherungs- 
weisen Integration der Differentialgleichung ebener, stationarer, laminarer 
Grenzschichten inkompressibler Fliissigkeiten bei vorgegebenen Randbedin- 
gungen entwickelt. Es diirfte den bisher angegebenen iiberlegen sein, da es die 
Vorteile eines geringen Arbeitsaufwandes und einer den Anforderungen ent- 
sprechenden Genauigkeit besser in sich vereinigt. Infolge seiner einfachen 
Handhabung ist es zur praktischen Verwendung geeignet. 

Die Entwicklungen schliessen an das Differenzenverfahren von H. GORTLER 
[1}?) an. Die vorgenommenen Anderungen zielen auf eine genauere Behandlung 
des wandnahen Gebiets. Die hierbei aufgewandte gréssere Sorgfalt ist gerecht- 
fertigt, da die Wand des umstr6mten Kérpers im Rahmen des Fortsetzungs- 
prozesses eine Linie der Unbestimmtheit darstellt. Wahrend die finiten Erset- 
zungen von H. GORTLER eine Verschiebung der Unbestimmtheitslinie in das 
Innere des betrachteten Gebiets bewirken, wird bei der neuen Methode im 
wandnahen Gebiet eine geeignete Polynom-Approximation vorgenommen und 
dadurch die fiir die numerische Durchfiihrung des Verfahrens besonders in 
der Nahe der Abldsungsstelle ungiinstige Verschiebung der Linie der Unbe- 
stimmtheit vermieden. Die neuen Approximationen bedingen bei jedem Fort- 
setzungsschritt die Durchfiihrung eines Iterationsprozesses, der theoretisch 
sehr durchsichtig und praktisch tiberaus einfach zu handhaben ist. Durch die 
Anderungen wird der Grundcharakter des urspriinglichen Verfahrens nicht ge- 
andert. Die Erhdhung des Arbeitsaufwandes — um 10 min pro Fortsetzungs- 
schritt, waihrend dieser bisher 1h in Anspruch nahm — ist unwesentlich und — 
wird durch die erzielte GenauigkeitserhGhung mehr als wettgemacht. Das Ver- 
fahren wurde an mehreren Beispielen erprobt ; dabei wurden — auch in der Nahe 
der Ablésungsstelle — gute Erfahrungen gemacht. 

Wie bei den Verfahren [1] und [2], so zeigt auch die nach dem verbesserten 
Differenzenverfahren errechnete Naherungslésung eine charakteristische Streu- 
ungserscheinung; in [9] wurde diese Aufspaltung der Lésung untersucht und 
eine zu ihrer Behebung geeignete Glattungsvorschrift angegeben. 

Es sei noch bemerkt, dass das Verfahren von K. SCHRODER([2], welches 
einen grosseren Zeitaufwand erfordert, mit der gleichen Genauigkeit arbeitet 


1) Mathematisches Institut der Universitat, Abteilung Angewandte Mathematik. 
*) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis auf Seite 397. 
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wie dasjenige von H. GOrTLER [1]. Auch hieriiber hat der Verfasser in [9] 
berichtet. 


§ 2. Das Randwertproblem der Grenzschichttheorie 


Die Geschwindigkeitsverteilung einer ebenen, stationaren, laminaren und 
inkompressiblen Grenzschichtstrémung genitigt bekanntlich den Prandtlschen 
Grenzschichtgleichungen. Bezeichnen x und y die Bogenlangen der angestrém- 
ten Wand bzw. der Normalen zur Wand, wu und v die zugehérigen Geschwindig- 
keitskomponenten, # den Druck, @ die Dichte und » die kinematische Zahigkeit, 
so lauten diese Gleichungen 


Ou Ou lap 07u 


a ie Oo, Q de |” dye” (fa) 
Ou Ov 
Rin nee (1b) 


Der Druck #(x) kann aus der als bekannt anzusehenden reibungsfreien Aussen- 
str6mung ermittelt werden). Bezeichnet u (x) die Langsgeschwindigkeit am 
ausseren Rande der Grenzschicht, so gilt in grenzschichttheoretischer Naherung 


Lape duy 
omar Ode 


(Ic) 


In die Gleichung (1) fiihren wir unter Verwendung einer charakteristischen 
Lange L, einer charakteristischen Geschwindigkeit U und der Reynoldsschen 
Zahl Re = UL/y dimensionslose Koordinaten ein 


7 x , Re , Pp 
a jae yay , p = 0 U?? | 
we (2) 
; U i re , U 
l= rien 2 ly =P. | 


Da im folgenden nur diese Koordinaten benutzt werden, kénnen die Striche 
wieder weggelassen werden. Nach Elimination von v mittels (1b)?) und # mit- 
tels (1c) lautet (1a) 


y © 
Ou ou f[ Ou duy 0>u 3a 
ia ee on OER ee (3a) 
0 


Diese Gleichung lasst sich durch Differentiation nach y in eine reine Differen- 
tialgleichung tiberfiihren. Es ist aber fiir das folgende zweckmiassig, sie in dieser 
integrierten Form zugrunde zu legen. 


1) Praktisch geschieht dies meist durch Messung des Druckes. 
2) Hierbei wird angenommen, dass die Wand undurchlassig ist: v(x, 0) = 0. 
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Die Funktion u(x, y) geniigt gewissen Randbedingungen. Wegen des Haf- 
tens an der Wand ist 
u(x, 0)=0. (3b) 


Am Rande der Grenzschicht geht die Langsgeschwindigkeit stetig in die der 
Aussenstrémung tiber. Diese Bedingung approximieren!) wir durch 


lim u(x, y) = u9(x). (3c) 
y—> 00 
Zur Lésung der Gleichung (3a) fiir x = x) unter den Bedingungen (3b) und 
(3c) bei bekanntem (x) ist noch die Kenntnis einer « Anfangsbedingung » 


u(x, y) = u(y) (3d) 


erforderlich; dabei kann #(y) nicht vollstandig beliebig sein?). 

Besonderes Interesse fiir die Praxis hat die sogenannte Ablésungsstelle 
(x4, 0), die durch du/dy (x4, 0) = 0 definiert ist. Da es zweifelhaft erscheint, 
dass die durch (3) bestimmte Funktion u(x, y) die (laminare) Geschwindigkeits- 
verteilung auch fiir x > x4 in geniigendem Masse approximiert, kann man 
sich bei der Lésung von (3) auf das Intervall x) < x < x4 beschranken. Wir 
nehmen fiir das Folgende an, dass die Lésung u(x, y) des Randwertproblems (3) 
in diesem Intervall existiert ; sie ist dann nach H. GORTLER [4] auch eindeutig 
bestimmt. 

Fiir die Lésung w(x, y) von (3) lassen sich einige weitere Bedingungen an- 
geben, die fiir y = 0 aus (3a) bzw. aus denjenigen Gleichungen folgen, die durch 
fortlaufendes Differenzieren nach y aus (3a) hervorgehen. Mit 


ed aie Gp araeh (4) 
lauten diese sogenannten Grenzschichtbindungen 
g(x) = p' = —Uuy Uy, 6s(x) = 0, Gq(X) == 2; Gy, | 5 
g(x) = 24,4, a(x) = 2a, ay, | 


Zur Integration von (3) hat man Naherungsverfahren anzuwenden, Liegt 
ein vorderer Staupunkt vor — diesem wird man dann die Koordinate x = 0 zu- 
ordnen —, so kann die Lésung u(x, y) in einem Gebiet 0 < x < XS hi; Vu 


') Genauer lautet diese Ubergangsbedingung u(x, d) =U o(*%), wenn 6 die Grenzschichtbreite 
bedeutet; die Art des Verhaltens der Lésung von (8) fiir yoo rechtfertigt jedoch nachtraglich 
die Approximation dieser Ubergangsbedingung durch (8c). ’ 

*) Die genaue Angabe der von u(y) zu erfiillenden Bedingungen, wie tiberhaupt der Randbe- 
dingungen des Problems (3), ist noch nicht méglich, da erst ein einschlagiger Existenzsatz diese 


Frage zu klaren vermag. Sicher ist jedoch lim u(y) = Uo(%p) erforderlich, wobei der Grenziibergang 
v—>oo 
in physikalisch sinnvoller Weise erfolgen muss. 
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aus wenigen Gliedern der Blasiusschen Reihe in einfacher Weise ermittelt wer- 
den ; die Grésse von x, richtet sich dabei nach der geforderten Genauigkeit und 
den vorliegenden Vertafelungen?). Dagegen hat man zur Berechnung der Lé- 
sung in %<%< x4 und damit der Abldésungsstelle x = x, — wie bei Stré- 
mungen ohne Staupunkt iiberhaupt — andere Verfahren anzuwenden. In § 3 
soll zunachst eine allgemeine Methode zur Gewinnung einer Naherungslésung 
aufgezeigt werden, wie sie durch die Struktur des Randwertproblems (3) nahe- 
gelegt ist. Die weiteren Ausfiihrungen dienen dann der Entwicklung eines 
speziellen Naherungsverfahrens dieser Art. 


§ 3. Fortsetzungsverfahren und Entwicklung der Problemstellung 


Ist fiir einen bestimmten Wert x = x* das « Profil» u(x*, y) bekannt, so 
lasst sich 0u/dx (x*,y) aus (3a) berechnen; unter Beriicksichtigung von 
du/dy (x*, 0) = 0 erhalt man nach L. PRANDTL {3} 


Ou Ou 1 02u i 1 / 024 , 
reas af ae (Sor + Mom) dy + (Som + Mot) (6) 
0 


Da nach (3d) das Profil u(x, y) bekannt ist, lasst sich w(x,, y) fiir hinreichend 
kleines x, — x) naherungsweise berechnen gemass 


Ou 
U(x, V) = U(%q, v) + (% — Xp) ne (Xo, ¥)- 


Ist %9 << 4% <+++ < %_ < %nay <+°* < %q eine geeignete Einteilung des Inter- 
valls <x), %4>, so erhalt man eine Naherungslésung von (3) durch rekursive Aus- 
wertung der Gleichungen 


0 
(Xn 41> Y) = Wns) + (ns — Fn) Se (Fn) (= 0,1,...). 


Allgemein spricht man von Fortsetzungsverfahren, wenn eine Naherungs- 
lésung von (3) durch « Fortsetzen » des fiir x = %) vorgegebenen Profils u(x, Y) 
_ mit Hilfe der Differentialgleichung (3a) unter Wahrung der Randbedingungen 
— gewonnen wird. Fiir diese hat (3d) die Bedeutung einer Anfangsbedingung. 

Die Gerade y = 0 stellt fiir den Fortsetzungsausdruck (6) eine Linie der Un- 
bestimmtheit dar; nach (3b) und (5) sind namlich fiir y = 0 der zweite Sum- 
mand sowie der Integrand des ersten Summanden unbestimmte Ausdriicke der 
Form 0/0. Da man beim Berechnen von Quotienten kleiner Zahlen eine ungin- 
stige Fehlerfortpflanzung zu erwarten hat, muss man bei finiten Approxima- 
tionen im wandnahen Gebiet erhéhte Sorgfalt anwenden. 

1) Die Koeffizienten in dieser Reihenentwicklung nach Potenzen von ~% setzen sich linear aus 


gewissen universellen Funktionen zusammen. Vertafelungen dieser Funktionen findet man fiir den 
besonders wichtigen Fall der zur Anstrémungsrichtung symmetrischen Strémungen bei A. ULRICH (5). 
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Unter den nach dem Fortsetzungsprinzip arbeitenden Verfahren haben sich 
die Differenzenverfahren bisher am besten bewadhrt. Bei diesen wird die Diffe- 
rentialgleichung (3a) in den Punkten (x;, y,) eines rechtwinkligen Gitternetzes 
der (x, y)-Ebene durch finite Naherungsgleichungen ersetzt. Diese Gleichungen 
lost man unter den gegebenen Randbedingungen auf, indem man die Bestim- 
mung der Werte u(x;, y,) mit 7 = +1 auf die der Werte mit 7 = » zuriick- 
fiihrt. Deshalb liegt auch bei den Differenzenverfahren eine ungiinstige Fehler- 
fortpflanzung im wandnahen Gebiet vor. 

Da die « Unbestimmtheit » von (6) fiir y = 0 keine Singularitat des Problems 
(3) darstellt, ist es — etwa durch Reihenentwicklung von (6) an der Stelle 
y = 0—miéglich, einen iiberall « bestimmten » Fortsetzungsausdruck anzugeben. 
Infolge der komplizierten Struktur der Differentialgleichung (3a) hat man 
jedoch bisher auf diese oder auf eine ahnliche Weise kein Fortsetzungsverfahren 
mit befriedigender Genauigkeit angeben kénnen. 

Es wird sich zeigen, dass die Lage der Linie der Unbestimmtheit von den 
vorgenommenen finiten Ersetzungen abhangt; das Gebiet einer ungiinstigen 
Fehlerfortpflanzung kann also verschoben werden. Die bisher bekannten Ap- 
proximationen, bei denen die Linie der Unbestimmtheit ausserhalb des Gebiets 
y = 0 liegt, haben einen zu grossen systematischen Fehler'). Bei den im Hin- 
blick auf ein einfach zu handhabendes Verfahren von H. GORTLER eingefiihrten 
finiten Ersetzungen liegt die Linie der Unbestimmtheit im Innern y > 0 des 
betrachteten Gebiets. Das ist fiir die praktische Rechnung — besonders in der 
Nahe der Ablésungsstelle — tiberaus ungiinstig, da hierdurch die Fehlerfort- 
pflanzung noch starker wird. Es liegt deshalb nahe, die Approximationen so 
abzudndern, dass einerseits die Linie der Unbestimmtheit wieder nach y = 0 
zu liegen kommt, andererseits aber der systematische Fehler und der Arbeits- 
aufwand nicht vergréssert werden. Indem wir ferner die Approximation in Wand- 
nahe méglichst genau vornehmen, tragen wir der Unbestimmtheit des neuen 
finiten Fortsetzungsausdrucks an der Wand Rechnung. 


In §4 sollen die fiir das Folgende wichtigsten Symbole eingefiihrt und die — 


Haupteigenschaften des Verfahrens von H. GOrTLER dargestellt werden. 


§ 4. Das Differenzenverfahren von H. Gortler 
Zugrunde liegt ein Netz von Gitterpunkten (x;, y,) mit den Koordinaten 


tps ty bt hy a(t OSE tsk Qala De) | 


. ma, 
Vyas bE etre Oo eee 1>o) | 


1) Ein derartiger Fortsetzungsdruck wird zum Beispiel bei K. SCHRODER [2], Seite 456, ange- 
geben; um die Genauigkeit der so ermittelten Werte zu erhéhen, sucht man diese Werte iterativ 
zu verbessern. Dieser Iterationsprozess konvergiert aber nach (9] gegen die nach H. G6rTxer [1] er- 


rechneten Werte, so dass die Linie der Unbestimmtheit letztlich doch wieder in dem hierfiit un- 
gunstigen Bereich y > 0 liegt. 
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und es wird zur Abkiirzung 
Us, n= W(%;, Ve) (8) 


gesetzt. Wir benutzen im folgenden die finiten Ersetzungen 


Ou Aik 

Ox (x;, Va) a Dh mit A; p= U%a4.4— U; 1k» (9a) 
Ou Vik 

ae Ca BT (eel) mwamite V4 b= Ui — Una (9b) 
07u Vien : * 

Oy? (%;, Yu) i. 4 [2 (k = 2) mit V2, = a V, k+1 Vega: (9c) 


Diese Approximation der Differentialquotienten durch Quotienten tibersprin- 
gender Differenzen ist im allgemeinen von grésserer Genauigkeit als die durch 
die tiblichen Differenzenquotienten. Auch fiir 


Va A Oy? (44, V1) (10a) 


lasst sich in einfacher Weise ein Naherungswert angeben!). 
Im Anschluss an (9) definieren wir noch die Gréssen 


Ou 
ir ane (10b) 
V2,= 42 —— - (445.0) Met aE (x,, 0) = — (up Up) e- (10c) 
Die Auswertung des Integrals 
ee 
ta=— | oe (ee 9) dy (11) 


erfolgt mit Hilfe der Trapezregel (Schrittweite /) bei Benutzen von (9a): 


Ror 1 Aes =~ 
aan -(3 ia fe 1 ai (21) (12) 


Damit lisst sich die Differentialgleichung (3a) in jedem Gitterpunkt (x,, Vx) 
mit k > 1 approximativ durch eine Gleichung zwischen Werten u,,,, ersetzen, 
Ist : 

yeas Vs 22 (0). 

1) Vgl. H. GOrTLER [1], Seite 175. Fiir das neue Differenzenverfahren werden wir in (18b) 

einen besseren Wert herleiten. 
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so kann man die Naherungsgleichung nach A; , auflésen und erhalt 


h k-1 
4 h(uo Uo)4 sip P Vier at Vir & Air 


1 
2 (Win iy Fn) 


Bei festem i = 7) und bekannten Werten 1; ,(k = 1) lasst sich hiernach Ai,,k 
unter Benutzung von (3b) rekursiv fiir k = 1, 2, ... ermitteln. Sind neben den 
Werten ;,, noch die Gréssen ;,_1,, (k 2 1) bekannt?), so ergeben sich nach 
(9a) auch die Werte u;,,1,,% (k = 1). Jede derartige Auflésung von (13) bedeutet 
die Durchfiihrung eines Fortsetzungsschritts. 

Die nach (13) ermittelte Naherungslésung erfiillt auch die Randbedingung 
(3c). Man folgert namlich?) aus (13) 


lim A;,, = 2 A(uo);- (14) 


k—oo 


: (k= 1) (13) 


Avn= 


Es gilt also mit lim 4; _ 1, ,= (ug); und lim u,, , = (up); auch lim #; , 1. = (Mo)i41> 
k-—>oo k—>0o k—> oo 
falls im Rahmen der Rechengenauigkeit 2 h(u9),; == (Uo);41— (o)«—1 ist. 

Das Auftreten des Summanden —V/;,;,/4 im Nenner von (13) zeigt, dass die 
Linie der Unbestimmtheit durch die finiten Approximationen in das Innere des 
Gebiets y = 0 verschoben ist. Wahrend der Nenner — als stetige Funktion in y 
aufgefasst — seine Nullstelle fiir +; << x4 im allgemeinen im Intervall 0 < y <1 
hat®), gut fir #,= x, 

1 


bal BY Semanal Vir. 


Entsprechendes lasst sich tiber die Nullstelle des Zahlers aussagen. Die 
Linie der Unbestimmtheit liegt also fiir x < x4 in der Nachbarschaft der Ge- 


1) Als Anfangsbedingung verlangt dieses Naherungsverfahren also die Kenntnis zweier benach- 

barter Geschwindigkeitsprofile; nach den Ausftihrungen der §§ 2 und 3 war zu erwarten, dass nur 
eine solche Anfangsbedingung bendtigt wird. Wir werden hierauf in § 6 kurz zuriickkommen. 

2) Man hat zu zeigen, dass (8c) fir x = X;41 erflillt ist, falls dieses fiir x = x, der Fall ist. 

Nun folgt aus lim u, ;,= (ug), die Existenz einer ganzen Zahl M, derart, dass fiir alle m BE M;, im 

k—0oo 
Rahmen der Rechengenauigkeit die Gleichungen Us, m = (Mo), und somit Vj, 41 = 0 und Vz m+2—=0 
m+1 m+ : 
gelten. Andererseits ist fiir jedes solche m 3’ A, , endlich, also V;,m429 5 4,;,,=0. Damit gilt 
=1 ral 
Ai, M+-2 : ¢ 
3) Fir #;< x4 ist V;,9 > 0 und somit u,; 


Yr 
= 2h(uo),; fir m= M,, das heisst (14). 
,0 — Vi, 0/4 < 0; andererseits gilt 
1 Ou 
2 


H 
)— Lye id xH(a 5) >0. 


— - 7. P 
Ma Kia wey, 
Dagegen ist fiir x; = x4 nach (5): 


a a ; 1 
u(x%4,¥) = y* + ri y®, das heisst u4 4 — <n Va,1 = 0(0°). 
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raden y = /. Dieses erklart die Tatsache, dass die numerische Ermittlung der 
wandnahen Werte «;, und auch uw, unbefriedigender verlauft als die der 
andern; das macht sich besonders in der Nahe der Ablésungsstelle bemerkbar. 

Die Verschiebung der Linie der Unbestimmtheit wird hier durch die Aus- 
wertung (12) des Integrals (11) verursacht. Es soll deshalb in § 6 eine andere 
Approximation des Integrals vorgenommen werden, bei der das vermieden 
wird. Zuvor werden in § 5 geeignete Approximationen von V9 und V7, her- 
geleitet. Einer besonders genauen Bestimmung dieser beiden Gréssen kommt 
grosse Bedeutung zu, da einerseits VY; 9 im neuen Differenzenverfahren eine zen- 
trale Stellung einnehmen wird, andererseits sich Fehler bei der Approximation 
von V;7, infolge der ungiinstigen Fehlerfortpflanzung in Wandnihe stark aus- 
wirken. 


§ 5. Finite Ausdriicke zur Bestimmung von /,, und /77, 


Zur Ermittlung von Naherungswerten der Gréssen V;9 und V7, approxi- 
mieren wir die Geschwindigkeitsverteilung u(x,;, y) im wandnahen Gebiet 
durch ein Polynom i,(y), das der Randbedingung (3b) sowie den ersten 
beiden Grenzschichtbindungen (5) geniigt: 


, Min ee haan, 15 
U(y) = boy + Be pag oo (15) 
Die noch freien Koeffizienten by, b,, ..., 6, legen wir derart fest, dass die Werte 
u(y,) (y =1,..., ) die vorgegebenen Werte u;, (y =1,..., 2) — in einer noch 


-naher festzulegenden Weise — approximieren. Dann wird auch eine befriedigende 
Ubereinstimmung in den Ableitungen der Funktionen u(x;, y) und %,(y) fiir 
kleine Werte von y vorliegen, so dass wir setzen kénnen 


On, 
a 16 
Pie CH! oral: (16a) 
eed yee 16b 
Viep=4l Oy? pel ( ) 


Wir haben im folgenden k = 5 und »=5 gewahlt; es hat sich namlich 
gezeigt, dass bei der fiir die praktischen Genauigkeitsbediirfnisse ausreichenden 
Festlegung eines Profils durch rund 15 Punkte dann die genaueste Bestimmung 
von V; 4 und V7; méglich ist. : 

Fordert man, dass die die Koeffizienten 0), 6, und db; festlegenden Beziehun- 
gen zwischen den Werten “,(y,) und w;,, linear sind, so folgt aus (15) und (16) 


5 
Vg = Ag Ven + 3A, tir, (17a) 
5 r=1 . 


5 
Poa = By V0 ce UW; y- (17b) 
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Die Koeffizienten A, und B, (r = 0, 1, ..., 5) sind von der Schrittweite unab- 
hangig. Die Linearitaét von (17a) wird sich in § 6 bei der Entwicklung des ver- 
besserten Differenzenverfahrens als sehr vorteilhaft erweisen. 

Zur Festlegung der Koeffizienten A, und B, fordern wir einerseits, dass die 
nach (17) fiir Potenzen 1., 2., 4. und 5. Grades?) von y ermittelten Werte mit 
den aus (10) resultierenden iibereinstimmen. Andererseits soll — der ungiin- 
stigen Fehlerfortpflanzung im wandnahen Gebiet wegen — der Einfluss von 
u, , und w, yin (17) ausgeschaltet werden. Eine einfache Rechnung liefert dann?) 


Vio 


ll 


—0,47347 V2, + 1,79091 u, 5 — 1,13331 «, 4 + 0,23210 u,,, (18a) 


2 
Vin 


I 


0,68972 V2, — 1,43631 w; , + 153392 u; , — 0,36535 u,,. (18b) 


Diese Formeln haben sich bei unseren Anwendungen sehr bewahrt. 


§ 6. Entwicklung des verbesserten Differenzenverfahrens 


Am Schluss von § 4 ergab sich die Notwendigkeit einer Approximation von 
(11), die keine Verschiebung der Tinie der Unbestimmtheit in das Innere des 
betrachteten Gebiets nach sich zieht. Zu diesem Zweck ersetzen wir die Funk- 
tion “(x,;, vy) in Wandnahe durch ihr Taylor-Polynom 3. Grades 


2 b 
hd 9) = weet Mihara ds (19) 


Fiihrt man die Differentiation 0u/d0x — gemiass (9a) — und die anschliessende 


Integration gliedweise durch, so ergibt sich unter Benutzung der bei dem Fort- 
setzungsschritt bekannten Grdésse 


Lip as Vv WEY ya (Vz. 1, Of Vena (20) : 


0 2 ‘< 
Vir=—| 5 (a, 9) dy=—(V, 0) (21) 


Da die Approximation (19) nur in Wandniahe giiltig ist, benutzen wir diese 
Quadraturformel nur fiir k = 1 und k = 2; fiir k => 3 dagegen verwenden wir 
weiterhin (12). Ersetzt man die Differentialquotienten wieder durch die Quo- 
tienten tiberspringender Differenzen (9), so ergeben sich aus (3a) die finiten 


1) Die Potenzen 0. und 3. Grades werden wegen (3b) und (5) nicht beriicksichtigt. 
*) Koeffizientenvergleich bei den genannten Potenzen in der Reihenentwicklung von (17a) 


bzw. (17b) an der Stelle y = 0 fiihrt auf je vier lineare Gleichungen fiir die sechs Unbekannten 
A, bzw. B,. A4y=4Ag=0 bzw. By = B,=0 liefert dann (18). 
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Gleichungen 

Meese. [iy Vito +h Gym) + 22), @=1,2) (2 

SF ap \ 41,0 ih) M4, oMolet Fe Vijng> (k= 1,2) (22a) 
1 k- 

A,,= CEL AMES Fas 3, ds. + 4h (ug Uo); + aA (& = 3) (22b) 


Die in (13) auftretende Verkleinerung von Zahler und Nenner, die sich be- 
sonders bei der Berechnung von 4; und A, unangenehm auswirkte, ist in 
(22) fiir diese beiden Werte vermieden. Andererseits hangt A, ;, fiir = 3 von 
den Gréssen A; und A; ab; hierdurch besteht, wie die Durchrechnung der 
Anwendungsbeispiele gezeigt hat, eine ausreichende Bindung zwischen den nach 
verschiedenen Vorschriften berechneten Werten A, , mit k = 1, 2 und A, , mit 
ez. 3. 

Aus (22) lassen sich die Gréssen 4;,, und damit u;,,, , mit festem 7 berechnen, 
falls neben den Werten w,;_,,;, und u,,;, (jeweils k 21) auch die Grosse V,,, 9 
bekannt ist. Dieses ist aber zunachst nicht der Fall; nach (17a) ergibt sich 
V;41,9 erst aus den Werten ;,,,;, gemass 


5 
V.1,0= Ao Vist Di 4, U; 49° (23) 


Diese weitere Beziehung zwischen den unbekannten u,,, ;, und V;,,,y neben den 
Gleichungen (22) erméglicht nun, den durch den Fortsetzungsprozess festge- 
legten Wert V;,,;,9 1m voraus zu berechnen. 

Es bestehen zwischen den sechs Gréssen V;,1,9, 4;,1, -.., 4;,5 sechs lineare 
Beziehungen: Die Gleichungen (23) sowie (22) fir k =1,...,5. Ersetzen wir 
in (23) “;41,, durch A; ;, gemass (9a), so ergibt sich mit der Hilfsgrésse 


K;=A NUS Ou Var) a Vag (24 
5 

Vi 9 ee DA, As, (25) 
(fil 


Die Gleichungen (22) sind offenbar von der Form 
Aia= Pia Vis1,0+ Qe: (26) 


wobei die Gréssen P,, und Q;,, wie auch schon K;, durch die Daten des 
(c — 1)-ten und 7-ten Profils festgelegt sind. Durch Einsetzen von (26) in (25) 
und Verwenden der Abkiirzungen 


Se eyes R= K+) A, OF (27) 
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erhalt man fiir V;,, 9 die lineare Gleichung 


V. 


7+1,0 


= SVis10+ T- (28a) 


Aus dieser ergibt sich V;,1,9, falls S;+1 ist, zu 


Vis1,0= [as z (28b) 

Die Berechnung von P, , und Q;,, und damit die von S; und 7; ist sehr um- 

stindlich; deshalb soll V;,,,9 anderweitig, und zwar mit Hilfe des Iterations- 
prozesses 


A’, = PVD + Qe (29a) 
Vero = K; +34, A), (29b) 


ermittelt werden. Zur numerischen Durchfiihrung dieser Iterationen bendtigt 
man die Gréssen P;,, und Q,,, nicht, da (26) nur eine das Wesentliche beto- 
nende Schreibweise von (22) war. Um die Konvergenz von (29) bei beliebigem 
V0, 9 Zu untersuchen, setzen wir (29a) in (29b) ein und erhalten 


Visio = SV ot T,. (30) 


Fir Zahlenfolgen mit diesem rekursiven Bildungsgesetz gilt, falls S;+ 1 ist, 


s s T; 
Vis0= ag + St (Mi se — as). (31) 
Somit konvergiert die Folge V7"), , fiir |S;| << 1; der Grenzwert ist dann unab- 


hangig vom Ausgangswert V7, ,, und zwar gleich der Lésung (28b) von (28a). 


Umgekehrt ist auch Konvergenz von (30) nur gegen diesen Wert moéglich. Aus 
drei aufeinanderfolgenden Werten 


ret ¥ V2.9 und Aci 
ergibt sich nach (30) 


Vint 6 Veen 
Sis rok eo (32) 
i+1,0 i+1,0 


Konvergenz der Folge V7", , gegen die Lésung (28b) von (28a) ist also gesichert, 
falls die Differenz aufeinanderfolgender Werte dem Betrage nach abnimmt. 


(32) zeigt ferner, dass fiir S; > 0 die Folge (30) monoton ist. 
Vermutlich gilt fiir x»< x,< x, 


0 syle (33) 
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Dieses war in allen durchgerechneten Anwendungsbeispielen der Fall. Ein 
Beweis der Giiltigkeit von (33) konnte jedoch bisher nicht erbracht werden?). 
Im folgenden wird (33) als giiltig angenommen. 

Da die Folge (30) unabhangig vom Ausgangswert konvergiert, kénnen im 
Verlaufe des Iterationsprozesses zur. Konvergenzbeschleunigung die errech- 
neten Naherungswerte V°;1} abgedndert werden: 


(s+1)* _ p7(s+1) +1 5 
yas = Vane ane; VET} a VA (34a) 


Einen Hinweis tiber die Wahl von o; liefert die aus (28b) und (30) folgende Be- 


ziehung 
Vaso = VALE Pte WES8 — 129, 0) (34) 


a 
4+1,0° "2+1,0 Sy 
a 


o, ist also von s unabhangig wahlbar. Zwar ist S;/(1 — S,) nicht bekannt; wegen 


k6énnen wir aber o; den Konvergenzverhaltnissen des vorhergehenden Fort- 
setzungsschritts entsprechend wahlen. Mit einiger Ubung wird man den Grenz- 
wert nach zwei, spatestens drei Schritten innerhalb der Rechengenauigkeit er- 
mittelt haben. 

Aus den Gleichungen (22) und (18) lasst sich bei vorgegebenem u(x) die 
Geschwindigkeitsverteilung u(x, y) angenahert berechnen, falls zwei Anfangs- 
profile «_,, und %»,, vorgegeben sind. Nach den Ausfiihrungen der §§ 2 und 
3 war zu erwarten, dass nur ein solches benotigt wird. Die vorliegende Erhéhung 


1) Nur fiir hinreichend kleines / liess sich bisher eine Aussage machen. Dann kann (19) in 
y 


0<y<5/als giiltig angesehen und somit (22a) fir k=1,...,5 verwendet werden. Nach (22a) 
und (26) ist 
Vi, Rk? 
Aes carat ete 
2 M5, & 


5 ; 
Durch Entwicklung von S; = » A, P;,, nach Potenzen ven / ergibt sich fiir aj+ 0 


7=1 
1 ao 
a A) 
si Wee + 0(/7) 
fur a, = 0 
: S,=1 +00). 


Es zeigt sich jedoch, dass die wahren Werte von S; — besonders infolge der speziellen Wahl der 4, 
in (18a) — wesentlich kleiner sind, als es die ersten Reihenglieder 


Pwartpn lassen: somit treten auch in der Nahe der Abldsungsstelle keine Konvergenzschwierig- 
< ? 


keiten auf. 
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liegt nun nicht daran, dass iiberhaupt finite Naherungen vorgenommen wurden, 
sondern daran, wie dieses geschehen ist ; denn sie wird vermieden, wenn man den 
Differentialquotienten 1. Ordnung du/dx (x;,y,) durch den Differentialquo- 
tienten 1. Ordnung (2; 1, — ™;,,)/h und nicht gemass (9a) durch den Differen- 
zenquotienten 2. Ordnung (u; 41, . — Ui—1,x)/2 ” ersetzt. 

Die Approximation (9a) wurde aus Genauigkeitsgriinden vorgenommen ; 
sie birgt jedoch — das haben die praktischen Anwendungen sowohl unseres Dif- 
ferenzenverfahrens als auch derjenigen von H.GORTLER{[1] und K. SCHRODER 
[2] gezeigt — in sich die Gefahr des Auftretens einer speziellen Streuungser- 
scheinung: Es zeigt sich némlich, dass die beiden Folgen errechneter Werte 
U_1h, U1, Ug,z,--- UN Uo ~, Uo,~, U4,p--- bei festem # fiir sich einen glatten 
Verlauf haben, sich jedoch mit wachsendem erstem Index mehr und mehr von- 
einander entfernen. Diese « Aufspaltung der Lésung » lasst sich, wie wir in [9] 
gezeigt haben, durch die Fortpflanzung von Fehlern in Strémungsrichtung er- 
klaéren und durch eine geeignete Glattung beheben. Wir kommen hierauf in 
§ 7 zuriick. 


§ 7. Die Praxis des Differenzenverfahrens 


Zur praktischen Anwendung des verbesserten Differenzenverfahrens wird 
man zweckmiassigerweise noch die folgenden Abkiirzungen einfiihren: 


| 2 Mn (k =1, 2) 
Neve (35a) 
1, af 
|2(™ 4-4 Via); (k = 3) 
7 h ° ji 9 2 
Zin = 4 (yg Uo); + se Vin = +e Vi»—Vi: ede BSh 
k-1 
See (t =3) (35) 


r=1 


Dann haben die Gleichungen (22) die iibersichtliche Gestalt 


1 Vi, 

i= | i41,0 — 44,1) aoe tt ; f (36a) 
1 

12> N; ~ Vis 0 L;,») V; 2+ Z; 2}, (36b) 
1 i] 

Ai. = Nae {Si5 Vi, + Z, 4} (k = 3) (36c) 


Sind zwei Anfangsverteilungen u_;, und wu, vorgegeben, so lasst sich das 
System (36) nacheinander fiir 7 = 0, 1,... unter Benutzung der Randbedingung 
(3b) auflésen. Nach (14) geniigt die so berechnete Naherungslosung auch der 
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Randbedingung (3c). Jede Auflésung des Systems (36) stellt die Ausfihrung 
eines Fortsetzungsschritts dar. Die numerische Durchfiihrung einschliesslich 
der iterativen Bestimmung von /,,, ) erfolgt an Hand eines Rechenschemas 
(vgl. das Beispiel in Tabelle 1). 

Die ersten zwei Spalten enthalten die Werte der beiden Ausgangsprofile 
u;_1,, und u,;;,. Die nachsten fiinf pes dienen der Berechnung der in (36) 
eingehenden Hilfsgréssen V;,, V2x, Vi 2/4, Ni, und Z;,. Um diese Spalten aus- 
zufiillen, bildet man zunichst aus den vorgegebenen Werten wu; , die tiber- 
springenden Differenzen 1. und 2. Ordnung: V,, mit k 21 und V7; mit k= 2 
und die Gréssen 


2 
Vir 


—4 1? (uy u%);, (10c) 


V, = —0,47347 V2, + 1,79091 u, , — 1,13331 w, , + 0,23210 4,,, (18a) 


1, 5? 


2 
Vi 


I 


0,689 72 V2, — 1,43631 w, , + 1,53392 u, 4 — 0,36535 u, ,. (18b) 


Aus diesen ergeben sich /,,/4 (k = 1, k = 3) und gemass (35a) und (35b) die 
Werte N; , und Z;,. Ferner berechnet man nach (24) mit Hilfe von (18a) bzw. 
nach (21) die Gréssen K;, L;, und L; 2: 


K,= =V, eyes oa Veron Vas (24) 
Li, a V;-10 — 0,16667 Vea — Ve 1:0) (20) 
Lis ae. Vi46 = 0,530.59 Vaan <7 Vu) (20) 


und vermerkt die Ergebnisse am Kopf des Rechenblatts. 

Nach diesen Vorbereitungen erfolgt die iterative Ermitttung von V;,1,0- 
Hierfiir ist die 3. Spaltengruppe vorgesehen, deren Elemente nicht vollstandig 
den k-Werten zeilenweise entsprechen. Die Iterationsvorschrift (29) lautet in 


der Terminologie dieses Paragraphen 


1 Waa 
(S) ee re Gy cs ila 
A { (72,0 [a Za (37a) 
Ss 1 Ss) E / 
AS, = ay (WE 0 — Fuad Mi + Zieh. (37b) 
1 ; 
AS, = ay USER Man + Zsah (k = 3, 4,5) (37c) 
k-1 
put Sie = Df Ave (k = 3, 4,5) 


Ve = 1/9091 A‘), — 1,13331 Aj), + 0,23210 gies (37d) 
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Die Konvergenz lasst sich durch Abainderungen der errechneten Werte VEx% 
Vii = Wha + a VET 8 — Vaso) (34a) 


mit 0; ~ 0;_, beschleunigen. Ausgehend von einem gerigneten Naherungswert?) 
V0", » = 1,5457 ergibt sich nach (37) eine neue Naherung Vo, = 1,5501. Dabet 
lassen sich die Gréssen Al!) , 2 i in einem Arbeitsgang auf der Rechenmaschine 
gewinnen; dasselbe ist der Fall fiir Al”, (k = 3, 4,5), falls zuvor die Grésse 
Sh (k= 3,4, ?) ermittelt ist. Es werden somit zwei Spalten benétigt, in denen 
die Werte S) und A‘, Aufnahme finden. Da eine Anderung von J), 4 gegen- 
aber Ve 6 pane. haben wir den errechneten Wert VY, » gemass (34a) korri- 
giert. Aus 2)" , = 1,5516 ergibt sich nach (37) eine neue Naherung 


V3), 9 = 1,5514. 


Um das Rechenblatt besser auszunutzen, haben wir die Werte des zweiten 
aa unter denen des ersten angeordnet. Im betrachteten Beispiel 
wurde /,°), , nochmals abgedndert und ein weiterer Iterationsschritt durchge- 
fiihrt. Mit Vorteil wird man auch die Werte V5") und Visa, auf dem Rechen- 
blatt vermerken; das kann in der Spalte der Grésse S\°), geschehen, da diese 
fiir k = 0,1, 2 nicht definiert ist. 

Die 4. Spaltengruppe dient der eigentlichen Berechnung der Grdésse A; , 
nach (36). Wieder sind zwei Spalten erforderlich, je eine fiir die Gréssen S; , 
und A; ;,. Sie sind wie die der 3. Gruppe zeilenweise auszufiillen ; das soll heissen, 
dass man nach 4, ;,(k 22) zunachst S;,,,= S;,,+4;,, zu berechnen hat. 
Erst dann lasst sich A; ;,., ermitteln?). 

Die letzte Spalte des Schemas nimmt die Werte 4; 41,4 = “;-1,x + 4;,x auf. 
Dabei zeigt sich gegebenenfalls die Notwendigkeit, dem Schema nach unten 
eine Zeile anzufiigen; das ist dann erforderlich, wenn der Anschluss an die 
aussere Geschwindigkeitsverteilung (i); ,, noch nicht hergestellt ist. 

Jeweils nach etwa vier Schritten wird man eine Glattung der errechneten 
Profile durchfiihren, um keinen Genauigkeitsverlust durch die in § 6 erwahnte 
Aufspaltung der Lésung zu erleiden. Das kann nach [9] auf folgende Weise 
geschehen: Man bildet aus den errechneten Werten 4, , die Gréssen 


1 


OF, ae & (4;,+ A; .1,%) (38) 


und tragt sie bei festem k iiber 7 auf. Es zeigt sich dabei, dass die Aufspaltung 


1) Diesen kann man einschatzen oder wie im Beispiel der Tabelle 1 durch lineare Extrapolation 
aus den Werten V;_,,9 und //, ,o geWinnen. — Hier und im folgenden geben wir des leichteren Ver- 
standnisses wegen die Z ahlenwerte des in Tabelle 1 betrachteten Beispiels mit an. 


*) Die Werte der fiinf ersten Zeilen der 4. Gruppe kann man unmittelbar vom letzten Itera- 
tionsschritt ibernehmen. 


201i 
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weitgehend behoben ist. Noch vorhandene Streuungen der einzelnen Werte 
lassen sich nun einfach graphisch glatten. Die so ermittelten Werte?) u, , werden 
im allgemeinen auch beziiglich ihrer Abhangigkeit vom Index k keine Streu- 
ungen mehr haben. Das kann man durch Auftragen der Gréssen 


Gs p= Up prim— Ur 


bei festem 7 iiber k oder durch Bildung des Differenzenschemas dieser Werte 
in bezug auf den zweiten Index itiberpriifen. 


Die Wahl der Schrittweiten # und / hangt von der geforderten Genauigkeit 
ab. Da eine Fehlerabschatzung bisher nicht durchgefiihrt werden konnte, lassen 
sich nur schwer allgemeingiiltige Aussagen machen. Bei einer geforderten Genau- 
igkeit von 1% wahlt man / etwa so, dass 15 Punkte die Anfangsprofile uv_,,, und 
Uo, festlegen; um einen Rechenschritt durchzufiihren, bendtigt man dann zirka 
70 min. Zweckmassigerweise wird dieser Wert / im Laufe der Rechnung festge- 
halten. 

Dagegen muss man die Wahl von A der Anderung der Geschwindigkeitsver- 
teilung u(*, y) bei wachsendem ¥ anpassen; insbesondere hat man h in der Nahe 
der Ablésungsstelle kleiner als sonst zu wahlen. Eine zu grosse Schrittweite macht 
sich auch in Unregelmassigkeiten der Folgen ..., #;_1,,, “;,, mit festem & bemerk- 
bar; wahrend die Aufspaltung bei normaler Wahl von /# nach vorgenommener 
Glattung erst langsam wieder auftritt, bewirkt eine zu grosse Schrittweite, dass 
gleich der erste neu berechnete Wert uw; ,,,, mit dem bisherigen (geglatteten) Ver- 
lauf ..., %j_1,%, Uz,, nicht in Einklang ist. 

Bei Schrittweitenverkleinerung ist ein neues zweites Ausgangsprofil erforder- 
lich. Soll etwa an der Stelle ¥ = x, die Schrittweite h durch h, ersetzt werden, so 
kann durch quadratische Interpolation das benétigte Profil an der Stelle 


#=4;-h=*3 
aus den (geglatteten) Profilen fiir 
*=H4,—-2Zh=%,_2, *=%,-h=4%x,_, und *=42%; 
gewonnen werden. Fiir h,=h/2 ergibt sich auf diese Weise 
1 3 3 
“irik = — gy Mayet Mie + > MR: (39) 


dabei wurde in Anlehnung an (8) u(4j;=q, Ve) = Ujpoa, p gesetzt. 


§ 8. Beispiele 


Im folgenden sollen Ergebnisse mitgeteilt werden, die bei der praktischen Er- 
probung des verbesserten Differenzenverfahrens erzielt wurden. Wahrend das 
zweite und dritte Beispiel nach den in § 7 gegebenen Anweisungen durchgefiihrt 


1) (38) liefert fiir das zuletzt errechnete Profil keine Korrektur; das ist nicht itberraschend, da 
eine mechanische Glattungsvorschrift auch das weitere Verhalten einer Kurve berticksichtigen 
muss. Legt man auf grosse Genauigkeit Wert, so wird man diesen Verlust in Kauf nehmen, denn 
gerade das letzte Profil hat als Ausgangsverteilung fiir die weiteren Fortsetzungsschritte besondere 
Bedeutung. 
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wurden, fanden im ersten Beispiel noch andere Ausdriicke fiir V;,9 und V2, sowie 
eine andere Glattungsvorschrift Verwendung?). 


1. Beispiel: Grenzschicht lings einer ebenen Platte mit linearer dusserey Geschwin- 
dighettsverteilung uy(x) = By — B, x 


Bei geeigneter Wahl der charakteristischen Gréssen, ndmlich JE TBE 
U=B,, nimmt die aussere Geschwindigkeitsverteilung in den dimensionslosen 
Koordinaten (2) die Gestalt 

Uo(¥) =1—-x% (40) 


an. L. HowartuH [6] hat fiir die Geschwindigkeitsverteilung der zugehérigen 
Grenzschicht eine exakte Lésung in Form einer unendlichen Reihe angegeben und 
die zur Berechnung der ersten sechs Glieder notwendigen Koeffizientenfunktionen 
vertafelt. Die Konvergenzgiite der Reihe gestattete jedoch keine genauere Aus- 
sage als 0,119 < ¥,< 0,129; deshalb wurden bei weiteren Berechnungen von 
Geschwindigkeitsverteilung und Ablésungsstelle diese Werte nur bis zu einer 
Stelle %,<*, benutzt. L. HowarrtH [6] erhalt so #¥,= 0,120, indem er fiir 
% > 0,0875 eine dem Pohlhausen-Verfahren verwandte Naherungsmethode an- 
wendet. H. GORTLER[1] benutzt fiir > 0,050 sein in § 4 dargestelltes Verfahren 
und erhalt 7, = 0,1185. 

Die Rechnung nach unserem Differenzenverfahren begann bei ¥ = ¥)= 0,050. 
Es wurde /=J)/0,032 gewahlt; das entspricht der Festlegung eines Anfangsprofils 
durch 10 Punkte. Durch diese Wahl von / war es moglich, die Ausgangswerte 


a4 


Differenzenverfahren 

— = = = GURILER 

Scale, ee HOWARTH 
Os ~ 


Fig. 1 
Grenzschicht bei linearer ausserer Geschwindigkeitsverteilung; Kurven @,(x). 


1) Fir V;,9 wurde der finite Ausdruck 
V0 = — 0,47697 V29 — 0,84518 wu, 1 + 1,355 70 u;, 9 


. 


+ 0,69035 us 3 — 0,67785 u;, 4 + 0,15482 u;, 5 
benutzt; er ergibt sich aus dem Polynom (15) mit k =5 durch die Forderung, dass die aus #,(y) 
gebildeten Ausdriicke V; 9, V;2, und V2 die vorgegebenen Werte annehmen. 2, wurde aus 

3 Toes 
V24= ws Veo + a, Vir2 


ermittelt: diese Formel erhalt man durch quadratische Interpolation in den Differenzen 2. Grades 
unter Beriicksichtigung von (5). Die Glattungen wurden graphisch vorgenommen. 
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bequem den Howarthschen Vertafelungen zu entnehmen. Im Intervall 
0,050 < * < 0,110 wurde mit h=0,005 gerechnet, fiir > 0,110 mit A = 0,0025. 
Dabei ergab sich a,(0,125)< 0. Eine genauere Auswertung lieferte 74= 0,124; 
es sei aber nochmals betont, dass nicht die Formeln (18) verwendet wurden. Mit 
(18) hatte sich vermutlich ein etwas kleinerer Wert ergeben. In Figur 1 sind zum 
Vergleich die Kurven a,(#) aufgetragen, wie sie sich'aus unserem Differenzenver- 
fahren, jenem von H. GOrTLER bzw. nach der Howarthschen Rechnung ergaben. 


2. Beispiel: Dey Hiemenzsche Kreiszylinder 


Von K. HremMENz[7] wurde die Druckverteilung an einem in Wasser 
(vy = 0,01 cm?s~4) getauchten Zylinder (Durchmesser 2 R = 9,75 cm) gemessen. 
Die Anstrémgeschwindigkeit war u.=19,2cms 1. Werden mit L=1cm, 
U =7,151 cms! dimensionslose Koordinaten (2) eingefiihrt, so lasst sich die 
Druckverteilung fiir 0 < * <7 in befriedigendem Masse durch die dussere Ge- 
schwindigkeitsverteilung 


Uo(x) = x — 0,006 289 x? — 0,00004615 x5 (41) 


wiedergeben. Die Approximation (41) kann somit zur Berechnung der Ablésungs- 
stelle nach dem Differenzenverfahren benutzt werden, falls die Rechnung *4< 7 
ergibt. Die von HIEMENz beobachtete Ablésung lag zwischen 6,81 und 6,98. 

Das Differenzenverfahren setzte bei = 4,5 ein. Es wurde /=0,3 gewahlt; 
dieses entspricht der Verwendung von 15 Profilpunkten. Die Werte der Anfangs- 
profile wurden aus den ersten ftinf Gliedern der Blasius-Reihe fiir zur Anstré- 
mungsrichtung symmetrische Str6mungen berechnet; dabei wurden die Vertafe- 
lungen [5] verwendet. Es wurden zehn Schritte mit h=0,2, dann je zwei mit 
h=0,1 und h=0,05 durchgefiihrt. Die Naherungsrechnung lieferte Ablésung 
unmittelbar hinter der Stelle # = 6,80. H. GORTLER[1] erhielt ebenfalls den Wert 
%4= 6,80, wahrend K. SCHRODER[2] %4= 6,87 ermittelte. Die erzielten Resul- 
tate sind in den Figuren 2, 3 und 4 veranschaulicht. 


a, 


Differenzenverfahren 
— — — GURTLER 
——---=- SCHRODER 


iL \ it a 
45 50 55 60 65 3 


Fig. 2 
: Grenzschicht am Hiemenzschen Zylinder; Kurven a(x). 
(Die von K. ScHrROpER ermittelten Werte waren nur kurz vor der Ablésungsstelle bekannt.) 
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Fig. 3 
Grenzschicht am Hiemenzschen Zylinder; Geschwindigkeitsverteilung 
bei konstantem Wandabstand. 


1 i 
O75 700 “CiK) 
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Fig. 4 
Grenzschicht am Hiemenzschen Zylinder; Auftragung einiger Profiie. 
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3. Beispiel: Dey Kreiszylinder mit potentialtheoretischer Druckverteilung 
Die Potentialstrémung zur Anstromgeschwindigkeit u., um einen Kreiszylin- 


der (Radius R) hat langs der Wand die Geschwindigkeitsverteilung 


(a) — 2 sins, (41a) 


falls mit L=R/2, U=u. dimensionslose Koordinaten (2) eingefiihrt werden. 
Die Geschwindigkeitsverteilung der zugehérigen Grenzschicht lasst sich fir 
0<x*< x, bei geeignetem %) wieder mit Hilfe der Blasius-Reihe ausrechnen. 


g 


7,0) 
Differenzenverfahren 
— — — BASIUS-ULRICH 


= 
0 yer, 
50 Be BY 3b ge + 


Fig. 5 
Grenzschicht am Kreiszylinder bei potentialtheoretischer Druckverteilung; Kurven a,(x). 
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Fig. 6 


Grenzschicht am Kreiszylinder bei potentialtheoretischer Druckverteilung; 
Vergleich einiger Profile. 
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Man approximiert (42a) durch das Taylor-Polynom 9. Grades 


1 i 1 
EE aT RV setae A Sealer Tel Bo} ie a 
um die Vertafelungen von A. UtricH[5] voll ausnutzen zu kénnen, A. ULRICH [8] 
benutzte die Werte der abgebrochenen Blasius-Reihe bis zur Ablésungsstelle hin, 
Durch Auflosung der Gleichung 0u/dy (¥4, 0) = 0 mit Hilfe der Vertafelungen [5} 
erhielt er 74 = 3,84,. 

Die Rechnung mit unserem Differenzenverfahren setzte bei # = 3,0 ein. Es 
wurde /= 0,3 gewahlt; dieses entspricht wieder der Verwendung von 15 Profil- 
punkten. Es wurden fiinf Schritte mit h=0,1, drei Schritte mit h = 0,05 und 
ein Schritt mit h=0,025 ausgefiihrt. Dabei ergab sich a,(3,675) < 0; eine ge- 
nauere Auswertung lieferte 74 = 3,67. Die erzielten Resultate sind in den Figu- 
ren 5 und 6 aufgetragen und mit jenen verglichen, die A. ULRIcH aus der Blasius- 
Reihe erhielt. 


LITERATURVERZEICHNIS 


[1] H. GOrTLER, Ein Differenzenverfahren zur Berechnung laminarey Grenzschich- 

ten, Ing.-Arch. 76, 173 (1948). 

[2] K. SCHRODER, Verwendung der Differenzenrechnung zur Bevechnung der lami- 

naven Grenzschicht, Math. Nachr. 4, 439 (1951). 

[3] L. PRaNDTL, Zur Berechnung der Grenzschichten, Z. angew. Math. Mech. 78, 77 

(1938). 

[4] H. Gorter, Uber die Lésungen nichtlinearer partieller Differentialgleichungen 
vom Reibungsschichttypus, Z. angew. Math. Mech. 30, 265 (1950). 

[5] A. ULtricu, Die ebene laminare Reibungsschicht an einem Zylinder, Arch. 
Math. 2, 33 (1949). 

[6] L. Howarth, On the Solution of the Laminar Boundary Laver Equations, Proc. 
Roy. Soc. London [A] 764, 547 (1938). 

[7] K. HrEMENz, Die Grenzschicht an einem in den gleichformigen Flissigketts- 
strom eingetauchten gevaden Kreiszylindey, Ding]. Polytechn. J. 326, 32 (1911). 

[8] A. Utricu, Die laminare Reibungsschicht am Kreiszylinder (1943); unver6f- 
fentlicht. Die wesentlichsten Resultate sind bei H. SCHLICHTING, Gyrenz- 
schichttheorie (Braun, Karlsruhe 1951), S. 128ff., angegeben. 

[9] H. Wittine, Uber zwei Differenzenverfahren der Grenzschichttheorie, Arch. 


Math. 4, 247 (1953). 


Summary 


The method of differences developed by H. GOrTLER[1] for the calculation of 
laminar boundary layers of plane, steady and incompressible flow functions in 
the neighbourhood of the point of separation hold only with a reduced degree of 
exactness. The same holds for the method of K. SCHRODER([2}. This is due to the 
displacement of a line of indetermination of the expression for the. continuation 
as a result of the finite approximations used. This defect is overcome by means 
of a polynomial approximation of the velocity distribution in the neighbourhood 
of the wall, and in addition the sensitivity of the method with respect to propa- 
gation of errors is generally reduced. The relatively slight increase in labour is 
more than offset by the increased degree of exactness. The splitting phenomena 
inherent in the above-mentioned methods may be avoided by means of a suitable 
smoothing process. The method has been tried with success on several examples. 


(Eingegangen: 16. Februar 1953.) 
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On the Equation for a Damped Pendulum under Constant Torque 


By WatLace D. Hayes, London?) 


The equation for a damped pendulum under constant torque, the same as one 
which is of interest in the study of synchronous motor oscillations, may be ex- 
pressed in the notation of SEIFERT®) 


a?6 dé e 
ees i a 1 
qe + a + sin@= 6B, (1) 


where « and f are positive real constants. The problem of interest is to determine 
the conditions under which a solution may or may not exist such that the quantity 


dé 
dy hge® (2) 


is non-negative and is periodic of period 2 2 in @. It is known?) that for B>1 
such a solution does always exist, and that for 0< #6<1 there exists an «(£) 
such that the conditions «< % and « > %) imply the existence and nonexistence, 
respectively, of such a solution. 
It is the purpose of this note to obtain bounds on a. We replace f by 6) accord- 
ing to the relation 
p=sind., 0>0.>> (3) 


< 


and replace 0 by 6 + a — 6); this translation does not affect the mathematical — 


problem. The differential equation is changed to one in y(@), yielding 


dy , 
y a +ay-+ sin(@+ a— 6)) — sind, = 0. (4) 


The procedure is to pick an approximate solution ¥(@) which is zero at 0= 0, 2a 
and to investigate the quantity g(@) defined by 


“ 0 a dy 3 #7 dy 
Saxe g(0) = y ar Sd le = aa ee | 
dy (5) 
=) 0 + « ¥ — cos @y sin@ — sin@y (1 — cos@) . | 


If this quantity is zero over the range 0 to 2.” then ¥(6@) is the solution desired 
and « = a. If it is negative or zero over the whole range then for y(0) = 0, 


y(2 ma) — y(0) = y(2m)>0. 


1) Office of Naval Research, U.S. Embassy, London. The opinions contained in this paper are 
the private ones of the writer and do not necessarily reflect the views of the United States Navy. 


2) GEORGE SEIFERT, On the Existence of Certain Solutions of a Nonlinear Differential Equation, 
ZAMP 3, 468-471 (1952). 
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Since it is possible to find a positive y(0) such that 
y(2 2) — y(0) <0 
the solution for some intermediate value of (0) will be periodic. Thus, 
g(0) 5 0(0S505 22) implies «Sa, (6) 


with equality only if (6) is itself a solution. On the other hand, if g is positive 
or zero over the whole range then for y(2 z) = 0, 


y(2 2) — y(0) = —y(0) <0. 


Consideration of the relation connecting any two solutions 


2a d0 
ye — 3 = C ex (-/—) 7 
Vier) 2 Pp ieee, (7) 
shows that a periodic solution with y(2 z) > 0 is impossible. Thus, 
g(0) 20 (0560522) implies «2 a, (8) 


with equality only if y(@) is itself a solution. 
The approximate solution is taken to be of the form 


y= jz) + cos 4) — 2 (0)| sins, (9) 


where /(0) = 1 and f is antisymmetric about 6 = a. The function g may be split 
into symmetric and antisymmetric parts 


g(8) 7 85 = ba? (10a) 
: 6 _. w so 
&,= 2 a|/ + cos Oy [2 — f cos — if sin=| — 2 sin 05 sin > (10b) 
2 ra 
g=5-[04+") cosy + DP = 21). (10c) 


The tests on g may be carried out in the interval 0 to wif both g, + g, and g, — &% 


are tested. 
For a lower bound on « we choose 


2 Sees See 
sin 6, = 2% ae cos6,, «? = |/3 cos?6,+ 1 — 2cos@), (11a) 


f,(0) = cos >, (11) 
i _. @ 
= SIs sin? — 4sin5], (alle) 
OH: at 
a, = —— sind sin. (11d) 


Since g, is negative for 0 << the test is made on g, — g,. Noting that «j < sin 4, 
we obtain ‘ ; ; 
8s, — & = 3 Sin I sin [sino + 3 sin > _ 4] <0. OS. ¢-S a) (12) 


Hence the quantity «, defined by equation (11a) is a lower bound to %». 
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For an upper bound to %) we choose 


as Se. 
sin 0) = a, a + cos6,, a= 2 sin) (13a) 
i 00 
f2(6) = 1 — sin (0o<60<2) 
; | (13b) 
=sing-—1, <0<2m | 
_ 8 6 
8s, = Sin A [2 = 2 cos 9 = sin 233 sin--cos-5 | § (13 :c)) 
3 : ng 7) shail 6 
La, = > 24 2 cos-$ Zein 3 sin-5 cos + 2 sin?jcos 5] . (13d) 


Since g,, takes both signs in the interval 0<@<2z the test must be applied to 
both g, + g, and g, — g,. Noting that «3/2 < sin 6), we obtain 
Bsa t+ Ba, = 48 sin? cos) =" : (0S0S2n)_ (14a) 


QO 


IV 


8s. = Sa. 


2 


<= 


a3 [2(2 — cos 5) + aoe cos 5] (: — sin) 20. (00S a) (14b) 


Hence the quantity «, defined by equation (13a) is an upper bound to a. This 
bound is considerably better than that of SEIFERT). 

The slope of the curve % (09) at 6) = 0 may be easily calculated for comparison 
with the same slope for the bounding functions obtained 


1 

ag(0) = —, (15b) 
4 

a{(0) = 1. (15c) 

It should also be noted that « (2/2) is finite, as 
a( 5) me (16a) 
n = 
aa(5-) = y2 . (16b) 


The solution in the vicinity of the neutral point 0 = 0 is always unstable, so 
that such a point (0 = 0, + 22m in SEIFERT’s notation) cannot be a limiting 
point for large ¢. The solution in the vicinity of the neutral point 6 = 2 6, — zis 
always stable for 6)< 2/2, so that such a point (0 = 0,-+ 227 in SEIFERT’S 
notation) can always be a limiting point for large ¢ with «> 0, and will be such 
a limiting point with «> a. 


1) GeorGE SrIFERT, On the Existence of Certain Solutions of a Nonlinear Dijjerential Equation 
ZAMP 3, 468-471 (1952). 


ee 


Vol. IV, 1953 Kurze Mitteilungen — Brief Reports — Communications bréves 401 


“ The case B = 1 or 6, = 2/2 has been avoided up to here, as has also the case 
% = %. If 6) = 2/2 and «< x, then the periodic solutions of the type considered 
do clearly exist. If 6, = 2/2 and « > « or if 0< 9,< 2/2 and « = a, the periodic 
solutions for y(@) exist, but are zero at 6 = 0, and the time ¢ given by the integral 


p= [ (17) 


does not converge. Such “‘periodic”’ solutions have no meaning in terms of the 
original function (z) and cannot be said to exist in the sense of the original intent 
of the problem. 

The problem of interest should be restated: to determine the conditions under 
_which a solution of equation (1) may exist such that the quantity d6/dt is positive 
_ and is periodic of period 2 z in 9. Then for B>1 such a solution exists, and for 
_0< 6 <1 there exists an «,(f) such that «< a and « > a imply the existence 
_ and nonexistence, respectively, of such a solution. For Bf = 0 and « = a = 0 any 

solution with y(0) > 0 is of this periodic type. 


Zusammenfassung 


Fiir die Differentialgleichung 


d*h dé : : Lg 
IE p/P SH = = een 
Geet hag sin} — sinh, = 0 (0<%< 7) 
gibt es eine Funktion «,(#,) mit den folgenden Eigenschaften: Wenn «< a, 
-existiert eine Lésung, worin d#/dt immer positiv und in @ periodisch ist. Wenn 
&% = %, existiert eine solche Lésung nicht. 
Man beweist, dass 4, < %< 0%, wo 


sin By = a V/o2 + 4 cosd, , a? = //3 cos?#@, + 1 — 2cos#,, 
1 ——— . O 

; x aa ‘0. 

sin, == a 2 Vo3 + 4cos8,, w= 2 sin->~- 


Die Grenze «, ist kleiner als die von G. SEIFERT gegebene. 


(Received: May 25, 1953.) 


Specific Heat of Gases at the Critical Point 
By John. F. Lee’), Raleigh, North Carolina 


Sonic velocity measurements have been used for some time to determine the 
specific heats of gases at low frequencies following the suggestion of EINSTEIN (it) 
that the reaction rates of reactive gas mixtures could be found from similar mea- 
surements. However, the equation for sonic velocity, expressed in terms of the 
ratio of the specific heats and the isothermal bulk modulus, becomes indeter- 
minate when applied to the critical point. In this paper a determinate expression 
for the sonic velocity at the critical point is derived which permits a direct solu- 


1) Department of Mechanical Engineering, North Carolina State College. 
2) Numbers in brackets refer to the Bibliography, page 404. 
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tion for the constant volume specific heat as a function of the critical values of 
temperature, specific volume, sonic velocity and Joule-Thomson coefficient. The 
calculated values obtained from this expression are examined in comparison with 
existing data. Finally, the approach to the problem from the standpoint of clas- 
sical thermodynamics is appraised in the light of some findings from statistical 
thermodynamics. 


Derivation of Cy, at the Critical Point 


The sonic velocity at zero frequency is expressed as follows: 


at = o— ($F) = - 2 ($F), (1) 


m \ Ov GNU 


where a = sonic velocity, m = molecular weight, v = molar volume, and —v(0p/0v), 
and — v(0p/dv) 7 are the isentropic and isothermal bulk moduli respectively. It 
can be shown that for any substance, 


T (0p/0T)2 


If this expression is substituted in equation (1), the following equation is obtained 
for the sonic velocity. 


2 vt (0p\2 v2 (op 
* mC, (57),- m (se). (2) 
At the critical point (0f/dv) ;=0 and equation (2) reduces to the determinate form 
ry Oe ea 
a -. mC, (se). (3) 


where a, is the sonic velocity at the critical point. 
Now the Joule-Thomson coefficient is defined as 


= (SF) 
which, with the aid of the second law of thermodynamics, can be shown to be 
equivalent to 
1 / Ov 
h=s Lig —o]. 5 
Cy OL. p ( ) 


The wellknown general thermodynamic relations, 


_ ee _ (0p/0T), 
OT]p (0p/dv) 7 


and 


A Ov Op 
Carte Sia pt aa 
» Cott (er) ar), 


when substituted in equation (5) yield, after simplification, the following equi- 
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valent expression for the Joule-Thomson coefficient 
Op) | Op 
oe Eadie Slee 
Op \2 Op\ — 
wE 
(ar).— Solar), 
Since (0p/0v) 7 = 0 at the critical point, equation (6) reduces to 


alt 
Me (Op/OT) (7) 


where “, = the Joule-Thomson coefficient at the critical point. 
When equation (7) is substituted in equation (3), the following expression for 
the constant volume specific heat at the critical point is obtained. 


C= = vy. (8) 


M \ Ae Me 


For the purpose of comparison some calculated values of the constant volume 
specific heats are listed with the zero pressure (no molecular interaction in the 
gas) specific heats in the following table. Data for the calculations and for the 
zero pressure specific heats were obtained from several sources ([2] to [7]). 


Constant Volume Molar Specific Heats at the Critical Point and at Zero Pressure 


C, [from equation (8)] Cc 


~ co 
Gas cal/mole degree cal/mole degree 
N, 11.3 4.98 
CO; 17.4 6.83 
O, eo 5-11 


The calculated values appear to be reasonable when compared with the zero 
pressure values of the specific heats and compare favorably with the values 
obtained from extrapolated P-v-T data. The calculated values at the critical point 
may be used to permit interpolation for intermediate values between the critical 
point and those points for which experimental data are available. It is believed 
that this is a more satisfactory method than that of simply extrapolating experi- 
mental specific heats to the critical point. 


Appraisal of Calculated Values of Cy at the Critical Point 


It is desirable to compare the calculated and calorimetric values of the spe- 
cific heats at the critical point, but unfortunately published values of the specific 
heats at the critical point are scarce. However, a calorimetric specific heats of 
approximately 50 calories per mole degree has been reported for CO, [8]. The diffe- 
rence between the calorimetric and calculated specific heats is to be expected in 
view of the limitations implied by the use of equation (8). This difference, it is 
believed, is of the correct order of magnitude to be accounted for by several fac- 
tors not anticipated by the classical thermodynamic approach employed thus far 
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The sonic velocity employed in equation (8) should be the thermodynamic 
value; that is, complete equilibrium must be attained. There is every reason to 
believe that the measured value of the sonic velocity is not the thermodynamic 
value because of the large sound attenuation at the critical point due to (1) scat- 
tering of the sound waves, (2) the increased absorption of high frequency sound 
waves at the critical point, and (3) configurational relaxations. 

Several investigators ([9] to [12]) have produced convincing evidence of vibra- 
tional heat capacity lag due to dispersion and high absorption of sound waves when 
a gas is rapidly compressed. Other investigators ([13] to [15]) have established the 
existence of rotational heat capacity lag. It is highly probable that when a gas is 
rapidly condensed the extensive clustering which ensues is affected by the sound 
waves causing a redistribution of the molecules in the clusters and/or a redistri- 
bution of the clusters accompanied by a configurational relaxation process. It is 
estimated from the previously mentioned papers that a total relaxation time of 
the order of 10~® seconds is required. 

From the foregoing discussion it must be concluded that (0f/0v) ; and (0p/0T),, 
are dependent on time. Therefore the system behaves as though it were only 
slightly compressible until sufficient time has elapsed for equilibrium to be 
established. 
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Zusammenfassung 


Es wird eine Ableitung einer bestimmbaren Gleichung fiir die akustische Ge- 
schwindigkeit am kritischen Punkt eines Gases gegeben, die gestattet, einen 
Ausdruck fiir die spezifische Warme zu gewinnen. Das Ergebnis wird hierauf 
vom Standpunkt der statistischen Thermodynamik erértert. 


(Received: February 4, 1953.) 
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————— ee ee 
Varia — Miscellaneous — Divers 


International Union of Crystallography 


The Executive Committee (General Secretary: Dr. R. C. Evans, Crystallo- 
graphic Laboratory, Cavendish Laboratory, Cambridge, England) has accepted a 
kind invitation from the French Government to hold the Third General Assembly 
and International Congress in Paris from 21 to 28 July 1954. These dates have been 
chosen in consultation with the French National Committee for Crystallography and 
with the National Committees of all the Adhering Bodies. It is hoped that this 
early notice will make it possible for crystallographers throughout the world 
to arrange to attend. 


Eine neue Standardfarbenkarte 


Der Fachnormenausschuss Farbe, ein Ausschuss des Deutschen Normenaus- 
schusses, hat in Zusammenarbeit mit dem Farbforschungslaboratorium des Ma- 
terialprifungsamtes Berlin-Dahlem in mehrjahriger Arbeit die experimentellen 
und technischen Grundlagen zu einer neuen Standardfarbenkarte erarbeitet, die 
in Zukunft die einheitliche Bezugsbasis fiir Farbenangaben in der Normung sein 
soll. Sie ist auf dem Gedanken der psychologischen Gleichabstandigkeit und 
Gleichwertigkeit der ausgewahlten Farben aufgebaut. Uber die wissenschaftlichen 
Grundlagen, die Ausfiihrungsform und die Anwendungstechnik in allen Zweigen 
der Praxis unterrichten die Aufsatze, die dariiber in der Zeitschrift «Die Farbe» 
(Verlag fiir angewandte Wissenschaften, Wiesbaden) in einem besonderen Heft 
(3/6, Bd. 1, 1953) zusammengestellt sind. Da die Arbeiten an der neuen Farben- 
karte unter der Leitung des Herausgebers der Zeitschrift, Prof. Dr. M. RIcHTER, 
gestanden haben, bietet das vorstehend erwahnte Heft Informationen aus erster 


Quelle. 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Description of a Magnetic Drum Calculator. By the Staff of the Com- 
putation Laboratory (Harvard University Press, Cambridge, U.S.A., 1952). 
318 pp., 196 figs.; $8.00. 

Es muss als iiberaus verdienstvoll bezeichnet werden, wenn sich die Erbauer 
einer programmgesteuerten Rechenmaschine nach Vollendung ihres Projektes 
der grossen Miihe unterziehen, eine detaillierte Beschreibung des Gerates auszu- 
arbeiten und zu publizieren. Der Stab des Harvard University Computation 
Laboratory hat dies mit dem vorliegenden Band 25 der «Annals of the Computa- 
tion Laboratory» zum drittenmal auf sich genommen, indem er zur elektronischen 
Rechenmaschine «Mark III» eine ausfiihrliche Darlegung des technischen, logi- 
schen und mathematischen Aufbaues verdffentlicht. Das Buch umfasst die elek- 
tronischen Schaltungen, die konstruktiven Einzelheiten, die logische Anordnung 
(wobei auch die mit dem verwendeten Zahlsystem zusammenhangenden Eigen- 
heiten erlautert sind), das verwendete Verfahren zur Berechnung der elementaren 
Funktionen, eine Anleitung zur Herstellung von Rechenplanen mit Beispielen 
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sowie eine Gebrauchsanweisung fiir das Gerat, und steht damit in seiner Art 
einzig da. — Gewiss verkérpert die beschriebene Rechenmaschine beziiglich ihrer 
elektronischen Schaltungen heute nicht mehr den neuesten Stand der Technik. 
Sie muss jedoch hinsichtlich ihrer Anpassung an die praktischen Bedirfnisse des 
Beniitzers als mustergiiltig angesehen werden. Dies gilt besonders fiir die An- 
fertigung der Zahlenstreifen und Rechenplane, die reichhaltigen Vorrichtungen 
fiir Eingang und Ausgang sowie die Bedienungs- und Kontrollmoéglichkeiten. ; 
Erwahnung verdient schliesslich die sorgfaltige Ausstattung des Buches. Die 
komplizierten Schaltbilder sind so aufgeteilt, dass sie nicht tiberladen erscheinen ; 
ganzseitige Photographien unterstiitzen die Anschauung. ASP. Speisey 


Fouriersynthese von Kristallen. Von W. Nowacki (Verlag Birkhauser, 
Basel 1952). 248 S., 120 Abb.; sFr. 34.30. 


Die Lagebestimmung der Atome in den Kristallen ist bei einfach gebauten 
Kristallstrukturen relativ leicht méglich. Sobald jedoch viele und auch verschie- 
denartige Atome in einem Kristallgitter vorkommen, wird die Ermittlung ihrer 
Schwerpunktslagen schwierig. In solchen Fallen wird bei der Kristallstruktur- 
bestimmung zuerst eine Symmetrieabklarung (Raumgruppenbestimmung) durch- 
gefiihrt, und zwar mittels Rontgeninterferenzaufnahmen, wobei Lage und Aus- 
wahl der vorhandenen gegeniiber den an und fiir sich médglichen Reflexen die 
Symmetriegruppe und die Metrik der Elementarzelle bestimmen lassen. Hierauf 
wird die Intensitaét der Reflexe zur Lagebestimmung der Einzelatome innerhalb 
der Elementarzelle beniitzt. Da bei Roéntgen- (und Elektronenbeugungs-) Auf- 
nahmen praktisch einzig die Elektronen des Kristalls streuend auf den Einfalls- 
strahl wirken, kann aus der Reflexintensitat einzig die Elektronenverteilung bzw. 
-dichte in der Elementarzelle bestimmt werden. Diese lasst dann ihrerseits die 
Atomlage ermitteln. Die Beziehungen zwischen der Elektronendichte in der Ele- 
mentarzelle und den Reflexintensitaten lassen sich als unendliche Fourierreihe 
darstellen. Dabei kénnen die Intensitaten (proportional zu den Amplituden der 
Streuwellen im Kristall) als Fourierkoeffizienten eingesetzt werden; die Phasen 
der einzelnen Streuwellen — in zentrosymmetrischen Kristallen zum mindesten 
die Vorzeichen der Koeffizienten — sind jedoch unbekannt. Die Methoden, welche 
diese Unbestimmtheit zu eliminieren gestatten, und die Rechenverfahren, welche 
die umfangreichen Reihenrechnungen beschleunigen kénnen, sind in den letzten 
Jahren von zahlreichen Autoren erweitert worden. W. Nowackr hat nun das in 
einer weitschichtigen Literatur zerstreute Material gesammelt und zusammenge- 
stellt. In seinem Buch fasst er die Rechenmethoden theoretisch zusammen und 
erlautert sie an Hand vieler Beispiele. Auch die den speziellen Verhaltnissen der 
Strukturbestimmung angepassten Patterson- und Buerger-Synthesen sind ein- 
gehend behandelt. Sodann bespricht der Autor die mechanischen und elektro- 
nischen Rechenhilfsmittel, das Lochkartenverfahren und optische Methoden, 
welche alle dazu dienen, die umfangreiche Rechenarbeit abzukiirzen bzw. zu be- 
schleunigen. Allen Abschnitten sind sehr eingehende Literaturverzeichnisse voran- 
gestellt, die es jederzeit erlauben, die Originalarbeiten zu Rate zu ziehen. Ein 
Bezugsquellennachweis fiir Rechenhilfsmittel erganzt das Buch. 

Das Werk von Nowackt ist die erste eingehende Zusammenfassung der spe- 
ziell fiir die Kristallstrukturbestimmung zugeschnittenen Fouriersynthese. Es 
fiillt damit eine bisher in Fachkreisen stark empfundene Liicke aus. Die zahl- 
reichen guten Figuren und drucktechnisch einwandfreien Tabellen erginzen die 
Arbeit aufs beste. Fiir Studierende, seien es Chemiker, Mineralogen oder Physiker, 
die ohne sehr eingehende Kenntnisse der Kristallographie an das Buch heran- 
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treten, ware eine kristallographische Einfiihrung sehr von Nutzen gewesen. Bei 
der rein mathematisch gehaltenen Auseinandersetzung mit dem Stoff findet der 
Leser auch keinerlei Beziehung zu dem in den iiblichen mineralogischen Lehr- 
biichern gebotenen Stoff. Es ist leider nicht dargelegt, wie man, ausgehend von 
einer Réntgenaufnahme, praktisch zu einer Fourierreihe gelangt, indem die 
Intensitatsbestimmung und -korrektur, die mit einer Reihe von kritisch zu bewer- 
tenden Messungen und Umrechnungen verbunden ist, gar nicht erwahnt wird. 
Der Zusammenhang des Experimentes — der Réntgenaufnahme — mit der mathe- 
matischen Behandlung ist damit fiir den nicht zum vorneherein Eingeweihten 
nicht gegeben. Stiinde eingangs ein entsprechendes Kapitel, so wire das vorlie- 
gende Buch auch eine sehr wertvolle Erganzung zu dem in der gleichen Reihe 
erschienenen Werk von E. BRANDENBERGER, Rdntgenographisch-analytische 
Chemie). Trotz dieser Mangel darf jedoch Nowackis Arbeit jedem Kristallogra- 
phen und Strukturchemiker aufs warmste empfohlen werden. W. Epprecht 


Magnetische Messungen an ferromagnetischen Stoffen. Von W. 
JELLINGHAUS (W. de Gruyter & Co., Berlin 1952). 163 S.; DM 18. 


Wenn ein Starkstromingenieur oder Elektrotechniker vor das Problem magne- 
tischer Messungen an ferromagnetischen Werkstoffen gestellt wird, ist sein Griff 
nach diesem einfiihrenden, klar und ausfiihrlich geschriebenen Buch als iiberaus 
gliicklich zu werten. Er findet im ersten Drittel des Bandchens einige einfiihrende 
Kapitel tiber MaB8systeme, magnetische Grundbegriffe, Erzeugung magnetischer 
Felder und Feldstarkemessung. Im zweiten Drittel ist die Messung der magneti- 
schen Induktion und der Magnetisierung in magnetischen Gleichfeldern beschrie- 
ben, und das letzte Drittel enthalt die Messmethoden mit niederfrequenten magne- 
tischen Wechselfeldern. Die Wirkungsweise und Handhabung der industriell ge- 
brauchten Eisenmessgerate ist tiberall sorgfaltig dargestellt. 

Hingegen hat der Autor bewusst die fiir die Schwachstrom- und Hochfre- 
quenztechnik ausserordentlich wichtigen Methoden der Permeabilitats- und Ver- 
lustmessung im Ton- und Radiofrequenzgebiet nur andeutungsweise eroOrtert. 
Die Behandlung der physikalisch aufschlussreichen Untersuchung der ferromagne- 
tischen Resonanzeffekte sowie der gyromagnetischen Effekte fehlt vollstandig. 
Diese Wiinsche der Fernmeldetechniker und Physiker ebenfalls zu erfiillen, ist 
freilich bei dem Umfang dieses Buches unméglich. Nichtsdestoweniger ist erfreu- 
lich, dass ein grosses und wichtiges Teilgebiet der ferromagnetischen Messtechnik 
eine leichtverstandliche, korrekte Darstellung gefunden hat, und es steht zu 
hoffen, dass dereinst auch die hier unbehandelten Untersuchungsmethoden 
ebenso gliicklich zusammengefasst werden. Hi. Labhart 


Finite Deformation of an Elastic Solid. By F. D. Murnacuan (John 
Wiley & Sons, Inc., New York, 1951). 140 pp.; $4.00. 

Das Buch, das in der von I. S. SoKOLNIKOFF herausgegebenen «Applied 
Mathematics Series» erschienen ist, behandelt die Elastizitatstheorie fester Kor- 
per unter Beriicksichtigung der Glieder héherer Ordnung in den Verzerrungs- 
grossen. 

Da sich der Verfasser konsequent der fiir diesen Zweck hervorragend geeig- 
neten Matrizenrechnung bedient, entwickelt er in einem ersten Abschnitt die 
spater bendtigten Grundformeln dieses Kalkils. . 

Der Aufbau der Theorie beginnt mit einer eingehenden Analyse der endlichen 
Verzerrungen, des zugehorigen Tensors und seiner Invarianten. Sodann wird der 
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7 


i 


auf die verzerrten Flachenelemente bezogene Spannungstensor eingefiihrt und — 


aus dem Prinzip der virtuellen Arbeit am deformierten Korper die allgemeine 


(auch fiir nichtisotrope Medien giiltige) Spannungs-Dehnungs-Beziehung her- 


geleitet, wobei die Massendichte der Formanderungsenergie als eindeutige Funk- 
tion des Verzerrungstensors angesetzt wird. 

Die Entwicklung dieser Funktion nach Potenzen der Verzerrungsgrossen bis 
zur dritten Ordnung fiihrt fiir den Spezialfall des isotropen Mediums, das zuerst 


untersucht wird, auf fiinf Elastizitatskonstanten, die ihrerseits noch von einem — 
hydrostatischen Ausgangs-Spannungszustand abhaéngen k6nnen. Hierbei zeigt 


es sich, dass im Bereich der endlichen Verzerrungen ein isotropes Medium, das 


anfanglich unter einem nichthydrostatischen Spannungszustand steht, bei einer 


weiteren Beanspruchung seine Isotropieeigenschaft verlieren muss. Die Theorie, 
die nun speziell auf die Verzerrung eines isotropen K6rpers unter hydrostatischem 


Druck angewandt wird, liefert eine im Rahmen der experimentellen Genauigkeit — 


vollkommene Bestatigung der Versuche, die 1948 von BRIDGMAN an Natrium 


unter Driicken bis zu 100000 Atmospharen ausgefiihrt worden sind. 
Die nichtisotropen Medien werden in einem besonderen Abschnitt behandelt, 
wobei sich der Verfasser im wesentlichen darauf beschrankt, die ersten drei Glie- 


der fiir die Reihenentwicklung der Formanderungsenergie bei verschiedenen 
: 


Formen kristalliner Medien zu ermitteln. 
Die beiden letzten Abschnitte sind weiteren Anwendungen gewidmet, und 


zwar der Scherung (wobei die auf den Ausgangszustand bezogenen Spannungen > 


eingefitihrt werden), dem einfachen Zug, der unter Aussen- und Innendruck ste- 
henden Kugelschale, dem entsprechend belasteten Rohr und der Torsion des 
Kreiszylinders. 


Das zum Teil aus Vorlesungen entstandene Werk ist in Form eines Lehrbuches — 
(mit zahlreichen Ubungsaufgaben) geschrieben und straff gegliedert. Es fordert 
vom Leser strenge Mitarbeit, zumal da ihm die Beweise fiir zahlreiche Einzelheiten — 


selbst tiberlassen werden. 
Die vom Verfasser angewandte Schreibweise der Matrizen kann bisweilen 


Missverstandnisse verursachen; auch ware es wiinschenswert, den Tensorbegriff— 


im Rahmen der Matrizenrechnung klarer herauszuarbeiten. Zu bedauern ist das 
fast vollige Fehlen von Literaturangaben. 


Nach Ansicht des Referenten erfiillt das Buch drei wichtige Aufgaben in 


hervorragender Weise, namlich erstens, dass es ein klares Urteil iiber Tragweite 


und Grenzen der klassischen Elastizitatstheorie erméglicht, zweitens, dass es 


einen wesentlichen Beitrag zu einer einheitlichen Auffassung dariiber liefert, wie 
die klassische Theorie auf eine solche endlicher Verzerrungen zu erweitern ist, 
und drittens, dass es dazu anregt, diese Theorie endlicher Verzerrungen nunmehr 
auch experimentell zu unterbauen. H, Kauderer 


Lezioni sulla teoria moderna dell’integrazione. Di M. Piconr e T. VIOLA 
(Edizioni scientifiche Einaudi, Torino, 1952) 403 p.; Lit. 5000.—. 


Les auteurs se sont proposé de donner dans cet ouvrage une exposition systé-— 


matique de la théorie moderne de l’intégrale en 1’édifiant sur les notions de fonc- 
tion additive d’intervalle et de fonction quasi-continue dans l’espace euclidien 
a un nombre fini de dimensions. L’ouvrage s’adresse donc en premier lieu aux 
mathématiciens qui désirent approfondir leurs connaissances sur Vintégrale de 
STIELTJES-LEBESGUE. Bien que sa lecture ne suppose que quelques connaissances 
préliminaires (résumées au chapitre I) sur les ensembles de points et sur les 
ensembles «ordonnés» d’opérations et qu’en général la démonstration des théo- 


: 
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rémes soit faite en détail, l’ouvrage sera surtout apprécié par un lecteur qui con- 
nait déja la théorie classique de LEBESGUE et est, par suite, 4a méme de constater 
Poriginalité de l’exposition en le comparant aux traités existants sur cette ma- 
tiére. Ci-dessous le résumé de son contenu: 

Table des matiéres. Introduction (p. 1 a 26). I. Ensembles de points et en- 
sembles ordonnés d’opérations (p. 27 a 46). II. Variation des fonctions d’intervalle 
(p. 47 a 66). III. Intégrale de RIEMANN-STIELTJES (p. 67 a 93). IV. Concept plus 
général d’intégration (p. 94 a 100). V. Masse lebesguienne des ensembles de 
points. Ensembles lebesguiens (p. 101 & 123). VI. Fonctions quasi-continues 
(p. 124 a 151). VII. Intégrale des fonctions quasi-continues et sommables (@. USZ 
a 181). VIII. Suites de fonctions sommables (p. 182 & 203). IX. Fonctions com- 
plexes. Convergence en moyenne. Applications & l’espace de HirBert (p. 204 
244). X. Réduction des intégrales multiples (p. 245 a 274). XI. Intégrale de 
‘LEBESGUE-STIELTJES relativement & une fonction de point (p. 275 a 314). 
XII. Fonctions additives d’ensembie (p. 315 a 340). XIII. Décomposition cano- 
nique des fonctions additives (p. 341 a 367). XIV. La dérivation des fonctions 
additives (p. 308 a 396). Index des notations, des exemples, des auteurs; indexe 
analytique (p. 397 a 402). M. Plancherel 


Hydrodynamische Grundlagen zur Berechnung der Schiffsschrauben. 
Von M. STRSCHELETZKY (Verlag G. Braun, Karlsruhe 1950). 257 S., 57 Abb.; 
DM 12.-. 

Die Schiffsschraube ist in ihrem tiblichen Aufbau mit wenigen, relativ breiten 
Fligeln ein hydrodynamisch sehr viel schwierigeres Objekt als etwa der Luftpro- 
peller. Wahrend man im letzteren Falle eine im ganzen befriedigende Theorie hat, 
ist dies bei der Schiffsschraube noch nicht der Fall. Der Verfasser hat sich nun 
bemiht, die komplizierten Vorgange der Rechnung zuganglich zu machen, und 
zwar geht er so vor, dass er das System der gebundenen und freien Wirbel passend 
approximiert und nun im Prinzip einfach, in der Durchfitithrung recht mtihsam 
nach Bior-SAavarRT die St6rungsgeschwindigkeiten ausrechnet. Leider konnte er 
die daraus entstandenen Zahlentafeln nur zum kleineren Teil dem Buche bei- 
legen. Es ist klar, dass auch dann noch sehr vieles durch Zwischenannahmen er- 
ganzt werden muss; aber es ist befriedigend, zu sehen, dass er diese moglichst ver- 
niinftig zu wahlen versucht. Insbesondere legt er grosses Gewicht auf Annaherung 
an die Minimalverteilungen. Die Fliigel werden als gebundene Wirbelflachen ein- 
gefiihrt, ihre Dicke durch passende Quellen- und Senkensysteme. Es ware zu 
wiinschen, dass die Ergebnisse einer solchen Rechnung mit Druckverteilungs- 
messungen an derart entworfenen Schrauben verglichen werden k6nnten. Erst 
dann kénnte man sich ein Urteil tiber den praktischen Wert der Theorie bilden. 
Auf alle Falle handelt es sich um einen ernst zu nehmenden Versuch, und man 
muss dem Verfasser fiir die ausserordentliche Miihe, die er sich genommen hat, 
volle Anerkennung zollen. J. Ackeret 


Grundlagen der Aeromechanik und Flugmechanik. Von A. PROLL. 
(Springer-Verlag, Wien 1951). 632 S., 278 Abb.; sFr. 49.50. 

Der Verfasser, der lange Zeit Professor in Hannover war, hatte kurz nach 
dem Ersten Weltkrieg eine «Flugtechnik» herausgegeben, die sich als ein recht 
niitzliches Buch erwies. Heute legt er uns ein ziemlich umfangreiches Werk vor, 
das ungefahr dasselbe Stoffgebiet behandelt. —- Was aber ist alles in den dreissig 
Jahren zwischen diesen beiden Ausgaben geschehen! Man versteht, dass auch auf 
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600 Seiten der ungeheuer angewachsene Stoff nur bruchstiickweise behandelt — 
werden kann. Es war zweifellos richtig, sich nunmehr nur auf die Grundlagen zu _ 
beschrinken. Aber selbst dies ist nicht leicht, da heute eine Armee von Forschern ~ 
und Ingenieuren an der Entwicklung des Flugzeuges arbeitet und man jetzt) 
Dinge zu den Grundlagen rechnen muss, die vor 10 oder 15 Jahren beinahe 
Terra incognita waren. Trotz redlichstem Bemiihen des Verfassers liegt uber 
dem Buch der matte Schimmer «der guten alten Zeit», in der PROLL selbst ja sehr — 
aktiv war. Manches, was heute in den Flugzeugfabriken zur Routine gehort, ist — 
nur knapp angedeutet, oder es wird von vorneherein auf die Literatur verwiesen — 
(allgemeine Profile, Fliigelschwingungen, neue Grundrissformen, Grenzschicht-§ 
umschlag, Rauhigkeitseinfltisse usw.). 

Das Buch gliedert sich, wie im Titel angedeutet, in zwei ungefahr gleich © 
grosse Teile. Im ersten kommen die Grundlagen der Aerodynamik, konforme Ab- : 
bildung, Profiltheorie und Grenzschichten und einiges aus der Gasdynamik zur © 
Sprache, im zweiten Teil in recht ausfiihrlicher, meist elementarer Form die Be- 7 
rechnung von Flugleistungen, Propellertheorie und Flugeigenschaften. Leider — 
beschrankt sich der Verfasser hier ganz auf die Flugzeuge mit kleinen Geschwin- © 
digkeiten und Motor-Propeller-Antrieb, so dass Mach-Einfliisse und die doch 
recht verschiedenen Eigenschaften des Strahlantriebes itiberhaupt nicht zur Be- : 
handlung kommen. 

Der Studierende wird gut tun, das Buch als eine niitzliche Zusammenfassung © 
der schon klassisch gewordenen Teile der Flugwissenschaft zu betrachten, wobei — 
ihm die sehr zahlreichen Beispiele zweifellos von Nutzen sein werden. Aber er 
wird nicht umhin kénnen, wesentliche Dinge aus anderen, heute ja ziemlich zahl- 
reichen Lehrbiichern diesen Grundlagen hinzuzufitigen. 

Ausstattung und Druck sind in gewohnter Weise vorziiglich. J. Ackevet 


«Die Farbe». Zeitschrift. (Verlag fiir Angewandte Wissenschaften GmbH., 
Wiesbaden). Einzelheft DM 7.80, Band (6 Hefte) DM 42.-. 


Die vierteljahrlich erscheinende Zeitschrift «Die Farbe», von der die ersten 
zwei Hefte vorliegen, ver6ffentlicht Originalarbeiten und zusammenfassende Ar- 
tikel, die sich mit der Farbe als optischer Erscheinung befassen. Als Herausgeber 
zeichnet Dr. MANFRED RICHTER, der durch seine zahlreichen Arbeiten auf dem 
Farbengebiet bestens bekannt ist und dessen Ubersicht iiber die Farbenlehre 
noch immer als deutschsprachiges Standardwerk auf diesem Gebiet gelten kann. 

Eine Arbeit von S. Réscu im ersten Heft der Zeitschrift befasst sich mit farb- 
metrischen Versuchen zur Papierchromatographie, wobei an Chromatogrammen 
einiger Aminosduren exakte Farbmessungen vorgenommen wurden und dabei vor 
allem der zeitliche Verlauf der Farbanderung registriert wurde. Ein Aufsatz von — 
A. KOHLRAUSCH liber das Arbeiten mit dem Helmholtz-Kénigschen Spektralfarben- 
Mischapparat wird vor allem die Physiker oder Physiologen interessieren, die 
selbst mit einem solchen Instrument arbeiten. 

F. Born referiert iiber die amerikanischen Normblatter iiber Farbmessung, 
wobei vor allem die Unterschiede zu den entsprechenden deutschen Normen dis- 
kutiert werden. In zwei weiteren Artikeln werden die Entwiirfe zu neuen Normen 
besprochen. Die erste Arbeit von W. ToELpTE befasst sich mit der Jodfarbskala, 
die durch Vergleich zur Kennzeichnung der Farbung von Lacken, Harzlésungen 
und dergleichen dienen kann, deren Farbe der einer Jod-Jodkalium-Lésung ahn- 
lich ist. Als Masszahl wird dann die Jodmenge in 100 ml einer waissrigen Jod- Jod- 
kalium-Lésung angegeben, die, unter genau vorgeschriebenen Bedingungen mit 
der zu untersuchenden Fliissigkeit verglichen, gleich hell erscheint. Von Vorteil 


a 
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ist dabei die genaue Reproduzierbarkeit der Vergleichslésungen; aber es kénnen 
natirlich nur Farbungen bewertet werden, die nicht allzustark vom Farbton der 
-Jodlésung abweichen. — Ein Artikel von F. Born bespricht den Entwurf zu einem 
neuen deutschen Normblatt iiber Farben und Farbgrenzen fiir optische Signale 
im Verkehr. ; 

Die Zeitschrift wird erganzt durch Buchbesprechungen und Nachrichten aus 
der Fachwelt. Hoffen wir, dass sie dazu beitrage, die moderne Farbenlehre, die 
ihre heutige festumrissene Gestalt erst in den letzten zwanzig Jahren gewonnen 
hat, immer mehr zum Allgemeingut zu machen. W. Grossmann 


Angewandte Radioaktivitét. Von K. E. ZImMEN (Springer-Verlag, Berlin 
1952). 124 S., 45 Abb.; DM. 18.80. 


Die Anwendung der natiirlichen und kiinstlichen Radioaktivitat hat in den 
letzten Jahren einen ausserordentlich starken Aufschwung genommen. Radio- 
aktive Isotope werden heute in der Medizin zur Bekampfung von Krankheiten 
sehr haufig angewandt, und die Zahl der Probleme, welche mit Leitisotopen in 
Chemie, Biologie, Physik, Technik usf. untersucht werden, ist in standigem 
Steigen begriffen. So ist es sehr zu begriissen, dass K. E..Z1MEN in gedrangter und 
fiir den Nichtphysiker doch sehr leichtfasslicher Art die Grundlagen der Radio- 
aktivitat sowie ihre Anwendung zusammengestellt hat. Mediziner, Techniker und 
Naturwissenschafter finden im ersten Teil eine iibersichtliche Erlauterung der 
Begriffe und Erscheinungen der Kernumwandlungen und -reaktionen und der 
dabei entstehenden Strahlungsarten. Im zweiten Teil ist eine grosse Zahl von 
Anwendungsbeispielen nach Fachgebieten gesondert zusammengestellt, so dass 
man eine Fille von Anregungen erhalt. Dieses Kapitel gibt auch dem Kernphy- 
siker einen guten Uberblick tiber die heutige Verwendung der Isotope. Viele 
Literaturzitate erméglichen es dem Interessenten, sich weiter in die theoretischen 
Arbeiten oder Anwendungsbeispiele zu vertiefen. Eine Reihe von Tabellen und 
Angaben tiber Strahlenschutz erganzen das Buch aufs wertvollste, so dass es 
jedermann empfohlen werden kann, der sich mit angewandter Radioaktivitat 
beschaftigt oder sich tiber deren Stand informieren mochte. W. Epprecht 


Theorie der Seometrischen Konstruktionen. Von L. BIEBERBACH (Verlag 
Birkhauser, Basel 1952). 162 S., 102 Abb.; sFr. 18.70. 


Die Theorie der geometrischen Konstruktionen befasst sich mit der Frage 
nach den Konstruktionsaufgaben, welche mit gegebenen Hilfsmitteln (zum Bei- 
spiel mit Zirkel und Lineal) losbar sind, oder etwa auch, welche Instrumente zur 
Lésung einer vorgegebenen Konstruktionsaufgabe erforderlich sind. Der sich 
scheinbar keiner klaren Methodik unterordnenden Eigenart des Stoffes entspre- 
chend, war bis heute eine eigentliche Theorie der geometrischen Konstruktionen 
ausgeblieben. Die in zahlreichen Lehrbiichern der Algebra eingeflochtenen Be- 
trachtungen iiber Konstruktionen beschranken sich in der Regel auf die Kon- 
struktionen mit Zirkel und Lineal. 

Im vorliegenden Buche fiihrt der Verfasser die typischen Fragestellungen der 
Theorie der geometrischen Konstruktionen in ihrer Wechselwirkung mit andern 
Zweigen der Mathematik vor. Insbesondere sind zahlreiche Briicken zur Algebra, 
zur Funktionentheorie und zur Zahlentheorie geschlagen. Das Hauptgewicht des 
Buches liegt weniger in der Behandlung von Einzelproblemen; das Ziel des Ver- 
fassers besteht vielmehr in einer Schilderung der Tragweite der verschiedenen 
Konstruktionshilfsmittel. Aus der Fiille der behandelten Konstruktionen mit 
vorgegebenen Hilfsmitteln seien genannt: Konstruktionen mit dem Lineal allein, 
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Konstruktionen mit dem Lineal allein nach Vorgabe eines gezeichneten regularen 
Polygons bzw. eines gezeichneten Kreises (Steinersche Konstruktionen), Kon- 
struktionen mit Zirkel und Lineal, Konstruktionen mit dem Einschiebelineal, 
Konstruktionen mit dem Rechtwinkelhaken und schliesslich Konstruktionen mit ~ 
dem Stechzirkel. 

Von den zahlreichen in die Theorie eingeflochtenen Einzelproblemen nehmen 
insbesondere die Winkeldreiteilung, die Konstruktion der regularen Vielecke und 
die Kreisbogenrektifikation breiteren Raum ein. Im Zusammenhang mit dem 
Fragenkreis der Winkeldreiteilung mit Zirkel und Lineal wartet der Verfasser mit 
einer Reihe von sehr netten und verfeinernden Resultaten auf. Die parallel mit 
der Bekanntgabe dieser neuen Ergebnisse laufende Polemik gegen einen bekannten 
Algebraiker wirkt allerdings nicht sehr verfeinernd. Nachdem BIEBERBACH — wie 
er selbst schreibt — mit seinem Buche am bewahrten Standpunkt festhalten will, 
dass jede Fragestellung und Meinung erlaubt sei, wirkt dieser Angriff etwas inkon- 
sequent. Der in Verbindung mit der Diskussion des Kreisbogenrektifikation — 
wiedergegebene funktionentheoretische Beweis der Transzendenz von e und x 
nach GELFOND und SIEGEL diirfte noch nicht sehr verbreitet sein. 

BIEBERBACHS temperamentvoll geschriebene Theorie der geometrischen Kon- 
stvuktionen vermittelt einen ausgezeichneten Einblick in die Vielfalt der Probleme 
und Methoden dieses Zweiges der Geometrie. Ganz besonders werden sich durch 
das vorliegende Werk die Lehrer der Mathematik angesprochen fiihlen, gibt es 
ihnen doch unzahlige Gelegenheiten zur Bereicherung und Belebung des Unter- 
richtes. Es bleibt nur zu wiinschen, dass auch die nie aussterbenden chronischen 
Winkeldreiteiler und Kreisquadrierer daraus ihren Nutzen ziehen kénnten. Das 
Buch ist als Band 13 der mathematischen Reihe der Lehrbiicher und Monogra- 
phien aus dem Gebiete der exakten Wissenschaften herausgekommen und ist vor- 
ziiglich ausgestattet. M. Jeger 


Advanced Statistical Methods in Biometric Research. By C. RapuHa- 
KRISHNA Rao (John Wiley & Sons, New York, 1952). 390 pp.; $7.50. 


Dieses Werk enthalt zunachst eine knappe Zusammenfassung der wichtigsten 
mathematischen Hilfsmittel: Vektoren, Matrizen, Determinanten, quadratische 
Formen. Hierauf werden die iiblichen Wahrscheinlichkeitsverteilungen angegeben 
und ihre hauptsachlichsten Eigenschaften hergeleitet. Der weitere Inhalt be- 
schlagt die folgenden drei Gebiete: 1. Probleme betreffend das Schatzen von 
Parametern und das Priifen von Hypothesen; 2. Probleme mit mehreren Varia- 
blen; 3. Probleme der Klassifikation. Zu jedem dieser Gebiete hat Rao in den 
letzten zehn Jahren wesentliche Beitrage geleistet. 

Das ausgezeichnete Buch von Rao kann jedermann empfohlen werden, der 
eine mathematisch strenge Begriindung der iiblichen statistischen Verfahren 
durcharbeiten und gleichzeitig einige der wichtigsten allerneuesten Entwicklungen 
kennenlernen will. Das Werk ist im Geiste von R. A. FISHER geschrieben: der 
Verfasser begniigt sich nicht damit, die mathematische Theorie zu erértern, er 
widmet auch der numerischen Auswertung besondere Sorgfalt, indem er ein 
insbesondere fiir mehrdimensionale Probleme ausserst wertvolles Rechenverfah- 
ren an verschiedenen Beispielen ausfiihrlich darstellt. 

Die von Rao gegebenen Methoden lassen sich nicht nur auf biometrische 
Fragen anwenden, wie man aus den von ihm gewahlten Beispielen und aus dem 
Titel des Buches schliessen kénnte; sie gelten ganz allgemein und werden zweifel- 
os in den verschiedensten Gebieten Anwendung finden, so insbesondere in der 
Industrie, woftir ibrigens schon jetzt in der Schweiz Belege vorliegen. A. Linder 
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«Das Buch zeichnet sich vor allem durch eine vor keiner Schwierigkeit zu- 
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